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0.1 Introduction 

Theoretical quantum optics studies "open systems," i.e. systems coupled to an 
"environment" ||l^, 0, pi H. In quantum optics this environment corresponds 
to the infinitely many modes of the electromagnetic field. The starting point 
of a description of quantum noise is a modeling in terms of quantum Markov 
processes [yj . From a physical point of view, a quantum Markovian description is 
an approximation where the environment is modeled as a heatbath with a short 
correlation time and weakly coupled to the system. In the system dynamics the 
coupling to a bath will introduce damping and fluctuations. The radiative modes 
of the heatbath also serve as input channels through which the system is driven, 
and as output channels which allow the continuous observation of the radiated 
fields. Examples of quantum optical systems are resonance fluorescence, where 
a radiatively damped atom is driven by laser light (input) and the properties 
of the emitted fluorescence light (output)are measured, and an optical cavity 
mode coupled to the outside radiation modes by a partially transmitting mirror. 

Historically, the flrst formulations were given in terms of quantum Langevin 
equations (as developed in the context of laser theory) which are Heisenberg 
equations for system operators with the reservoir eliminated in favor of damping 
terms and white noise operator forces (see references in [|]). The alternative 
and equivalent formulation in terms of a master equation for a reduced system 
density operator together with the quantum fluctuation regression theorem, has 
provided the most important practical tool in quantum optics, in particular for 
nonlinear systems 0, 0| . 

The rigorous mathematical basis for these methods is quantum stochastic 
calculus (QSC) as formulated by Hudson and Parthasarathy m. QSC is a non- 
commutative analogue of Ito's stochastic calculus |^. The basic ingredients 
are "white noise" Bose fields b{t), 6(t)^ with canonical commutation relations 
[b{t),b{t'y] = d{t — t'). In quantum optics these Bose fields can be consid- 
ered as an approximation to the electromagnetic field. The connection between 
these abstract mathematical developments and the physical principles and foun- 
dations of quantum optics is discussed in particular in the work of Barchielli 
[RJ, |, 111, and Gardiner and coworkers ||l|, |l|, H). Gardiner and CoUett 
[I2| showed the connection between the more physically motivated quantum 
Langevin equations and the more mathematically precise "quantum stochastic 
differential equations" (QSDE). Furthermore, these authors introduced as an 
essential element of the theory an interpretation of the Bose-fields as "input" 
and "output" fields corresponding to the field before and after the interaction 
with the system. QSC allows also the development of a consistent theory of pho- 
todetection (photon counting and heterodyne measurements) |0, 04] in direct 
relationship with the theory of continuous measurement of Srinivas and Davies 

While most of the theoretical work in quantum optics has emphasized quan- 
tum Langevin and master equations M, recent developments and applications 



have focused on formulations employing the quantum stochastic Schrodinger 
equation (QSSE) Jp JlC , [l^, pi| and its c-numbcr version, a stochastic Schrodinger 



equation (SSE) ]2|, 10| , 
29| 
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27,28, 
Ol |31|, |33, M. Phenomenological SSEs have been given in Refs. |34|, 35, 36, 
37 1 . The QSSE is a QSDE for the state vector of the combined system + heat- 
bath with Bose fields h{t) and 6(i)^ as noise operators. This equation generates 
a unitary time evolution, and implies the master equation and the quantum re- 
gression theorem for multitime averages, and the Quantum Langevin equation 
for system operators as exact results. 

The SSE, on the other hand, is a stochastic evolution equation for the sys- 
tem wave function with damping and c-number noise terms which implies a 
non-unitary time evolution. The relationship between the QSSE and the SSE is 
established by the continuous measurement formalism of Srinivas and Davies for 
counting processes (photon counting) |]1^ . This continuous measurement theory 
has the interpretation as a probabilistic description in terms of quantum jumps 
in a system. Thus the continuous measurement theory can be simulated proba- 
bilistically, and this simulation yields a sequence of system wave functions with 
jumps at times ti,t2, ■ ■ ■ , t„ by a rule which is directly related to the structure of 
the appropriate QSDE. The trajectories of counts generated in these simulations 
have the same statistics as the photon statistics derived within the standard pho- 
todetection theory, and in this sense the individual counting sequences generated 
in a single computer run can, with some caution, be interpreted as "what an 
observation would yield in a single run of the experiment." These simulations 
can thus illustrate the dynamics of single quantum systems, for example in the 
quantum jumps in ion traps |p8| or the squeezing dynamics M. In continuous 
measurement theory the time evolution of the system conditional on having ob- 
served (or simulated) a certain count sequence is called a posteriori dynamics. 
The master equation is recovered in this description as ensemble average over 
all counting trajectories. Again, in the language of continuous measurement 
theory, this corresponds to a situation where the measurements (the counts) are 
not read, i.e. no selection is made (a priori dynamics). Simulation of the SSE 
and averaging over the noise provides an new computational tool to generate 
solutions of the master equation. An important feature, first emphasized by 
Dalibard, Castin and M0lmer ^^ is that one only has to deal with a wave func- 
tion of dimension N , as opposed to working with the density matrix which has 
N'^ elements. Thus, simulations can provide solutions when a direct solution of 
the master equation is impractical because of the large dimension of the system 
space. 

When instead of direct photon counting we perforin a homodyne experiment, 
where the system output is mixed with a local oscillator, and a homodyne 
current is measured, the jumps are replaced by a diffusive evolution. These 
diffusive Schrodinger equations were first derived by Carmichael M, H from his 
analysis of homodyne detection, and independently in a more formal context by 
Barchielli and Belavkin |11| . Somewhat earlier equations of this form have been 



postulated in connection with dynamical theories of wave function reduction 

In these notes we will review in a pedagogical way some of the recent devel- 
opments of quantum noise methods, continuous measurement and the stochastic 
Schrodinger equation and some applications, and the main purpose is to em- 
phasize and illustrate the physical basis of such a formulation. 

0.2 An introductory example: Mollow's pure 
state analysis of resonant light scattering 

As an introduction and historical remark we briefly review Mollow's work on 
the pure state analysis of resonant light scattering [^ . Twenty years ago Mol- 
low developed a theory of resonance fluorescence of a strongly driven two-level 
atom, where he showed that the atomic density matrix could be decomposed 
into contributions from subensembles corresponding to a certain trajectory of 
photon emissions, which could be represented by atomic wave functions. This 
formulation anticipated and provided the basis O, O, pq] for some of the re- 
cent developments in quantum optics from the Stochastic Schrodinger Equation 
for system wave functions to our understanding of the photon statistics in the 
framework of continuous measurement theory. Below we will briefly review some 
of Mollow's ideas in a more modern language, closely related to the discussions 
of the following sections. 

Hierarchy of equations 

We consider a two-level system with ground state \g) and excited state |e) 
which is driven by a classical light field and coupled to a heatbath of radiation 
modes. The starting point of the derivation of Mollow's pure state analysis is 
the definition of a reduced atomic density operator in the subspace containing 
exactly n scattered photons according to 

p(«)(t) = TrB{p(")p(t)} (n = 0,l,2,...) (1) 

with pit) the density operator of the combined atom -f- field system, P^") the 
projection operator onto the n-photon subspace, and Tr^ indicating the trace 
over the radiation modes. The probability of finding exactly n photons in the 
field is 

P(") (i) = Tr5p(") (t) ^ pI^J (t) + (^ (t) (2) 

with Trg the trace over the atomic variables (the "system" ) . It has been pointed 
out in Refs. |^, Q Q that p^") (t) obeys the equation of motion 



Here 

H,B = {-A-t^r)a,e-(^na++H.c\ (4) 

is the non-hermitian Wigner-Weisskopf Hamiltonian of the radiatively damped 
and driven two-level atom with a+ — \e){g\, U- — |(?)(e|, az = |e)(e| — \g){g\ = 
o'ee — o'gg the Pauli spin matrices describing a two-level atom, A the laser 
detuning, F the spontaneous decay rate of the upper state of the two-level 
system and il the Rabi frequency. Summing over the n-photon contributions 

oo 
n=0 

with p{t) the reduced atomic density operator, Eq. (ph reduces to the familiar 
optical Bloch equations (OBEs) 



It 



p{t) = -z (iJeffP(i) - p{t)Hl^) + ra_p(t)a+ . (6) 



The n-photon density matrix p'"' [t) is seen to obey a hierarchy of equations 
where the n — 1 photon density matrix provides the feeding term for the n- 
photon term, . . . ^ p^'^^^'> -^ p^"^ .... By formal integration of this hierarchy 
we obtain Mollow 's pure state representation of the atomic density matrix 
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p(t) = |^c(i)>('/'cWI +1] / rf^n / dtn-1... dti (7) 

n=l"'0 Jo Jo 

IV'cCilin, t„_i, . . . , ti)) {^l}c{t\tn,tn^l, . . . , tl)| (8) 

with atomic wave functions obeying the equation of motion 

-\Mt\tn,-.-,h)) =-iH,s\M-t\tn,.-.,h)) (i>t„>t„-l >...>tl) (9) 

and defined recursively by the jump condition 

\lpcitn\tn,tn-l,-. .,tl)) = Vf Cr_ |'0c(i„ |i„-l , . . . , il)) (10) 

with \\ipc{t = 0)1 = 1. Eq. {uu gives the atomic density matrix in terms of pure 
atomic states ipdtltn, ■ ■ ■ , ti). 

Interpretation 

The essential features behind this construction are as follows. According to Mol- 
low [|9| Eq. ( 0) can be interpreted as the time evolution of an atom in a time 
interval (0, t] which emits exactly n photons at times ti, t2, . . . , t„. According to 
(lOh atoms in the ground state with n photons in the scattered field are created 
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Figure 1: Plot of a realization of the Monte Carlo wavefunction as a function of 
time: excited state probability |(e|'0c(i))P ]l9[ . 



at a rate rla-Vcl^Kn-i, ■ ■ ■ ,ii)|P- Each photon emission is accompanied by 
a reduction of the atomic wave function to the ground state \g) as described 
by the operator cr_ in Eq. (IIQ). This is what we call a quantum jump. The 
time evolution between the photon emissions is governed by the non-hermitian 
Wigner-Weisskopf Hamiltonian Hgg which describes the reexcitation of the ra- 
diatively damped atom by the laser field. 

As has been pointed out in RefJSSi in the context of discussion of quantum 
jumps in three-level systems |^, |4l|, Eq. (|^) has the structure expected from 
the Srinivas Davies theory of continuous measurement for counting processes 
K^ , and in particular supports the interpretation of 
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,h) = \\M-t\tr, 



(11) 



as an exclusive probability density that the atom emits n photons a times 
ii,i2, ■ ■ ■ ,tn (and no other photons) in the time interval (0,i]. Thus the time 
evolution can be simulated: Fig. Qshows the excited state population | (e|?Ac(0) P 
with normalized wavefunction ■(/'c(i) = V'c(i)/||V'c(i)||, corresponding to a single 
run as a function of time for a Rabi frequency f2 = F and detuning A — —T. 
The decay times (indicated by arrows) and quantum jumps where the atomic 
electron returns to the ground state are clearly visible as interruptions of the 
Rabi oscillations. After a quantum jump the atom reset is to the ground state, 
and it then is reexcited by the laser: in the photon statistics of the emitted light 
this leads to antibunching, i.e. two photons will not be emitted at the same time 
[^. Averaging over these trajectories gives the familiar transient solution of the 
optical Bloch equation. Although MoUow's derivation and somewhat intuitive 
interpretation was restricted to a specific example, the basic structure and ideas 
which emerge from his analysis are valid in a much more general context which 
is what we will be concerned with in the following sections. 



0.3 Wave-function quantum stochastic differen- 
tial equations 

0.3.1 The Model 

The standard model of quantum optics K, nfl] considers a system interacting 
with a heatbath consisting of many harmonic osciUators representing the elec- 
tromagnetic field. The Hamiltonian for the combined system is (?i = 1) 

-fftot = Hsys + Hb + Hint (12) 

with 

Hb = f dajLubiuj^biLj) (13) 

the Hamiltonian for the heatbath, and b{uj) the annihilation operator satisfying 
the canonical commutation relations (CCR) 

[b{ij),b{ijy]=Siu;-io') . (14) 

For simplicity, we assume only a single heatbath. Generalization to many reser- 
voirs represents no conceptual difficulty. The interaction Hamiltonian (nsf) is 
based on a linear system - field coupling in a rotating wave approximation 
(RWA), 

Fint = i-^ f duj K{uj)[cb{uj)^ - 6(w)ct] (15) 

with c a system operator (the "system dipole"), and k(w) coupling functions. 
The system Hamiltonian i/gys is left unspecified. 

Approximations 

Quantum noise theory is based on the following approximations ||l|: (i) Rotat- 
ing wave approximation (RWA) and smooth system-bath coupling, and (ii) a 
Markov (white noise) approximation. 

• Rotating wave approximation: For simplicity, we assume that for the bare 
(interaction free) system the system dipole c evolves as c(i) = ce~'"°* 
with LjQ the resonance frequency of the system. Thus the system will be 
dominantly coupled to a band of frequencies centered around ljq. Validity 
of the RWA requires that the frequency integration in iJjnt is restricted to 
a range of frequencies 

PUJ0+-8 

/ duj... (16) 

with cutoff d -^ LOf). This assumes a separation of time scales: the optical 
frequency ojq is much larger than the cutoff ^ which again is much larger 



than the typical frequencies of the system dynamics and frequency scale 
induced by the system - bath couplings (decay rates). Furthermore, we 
assume a smooth system-bath coupling k(uj) in the frequency interval 
(EgI): we set k{ll)) ^- 1 (a constant factor can always be reabsorbed in a 
definition of c) . 

We transform to the interaction picture with respect to the free dynamics 
of the "bare system + field." In the interaction Hamiltonian (p^ ) this 
amounts to the replacements c -^ ce~*'^"*, 5(a;) -^ 6(ti;)e^*"*, so that in 
the interaction picture the interaction Hamiltonian becomes 

^iLt^W=*rn0^c-6(-)(t)ct] (17) 

with 

b^^\t) = ^= / duj 6(w)e-'("-"°)*. (18) 

V27r Jcjo-1? 

The time evolution operator [/'^''^(t) from time to i in the interaction 
picture obeys the Schrodinger equation 

^^U^^\t)=-^{H + Hl^,\t))U^^\t) (19) 

with H the transformed system Hamiltonian. Note that by going to the 
interaction picture "fast" optical frequencies loq have disappeared ( "trans- 
formation to a rotating frame" ) . 

The Markov approximation or white noise approximation: this consists of 
taking the limit -i? ^ oo in Eq. (|l9|), 

b^^\t) ^ b{t) -.^ ^= / dt^6(w)e-'("-'^°)* . (20) 

so that the commutator for the fields b{t) in the time domain acquires a 
^-function form, 

[b{t),bity]=S{t-t'). (21) 



The operator b{t) is a driving field for the equation of motion (^^ at time 
t. In the following we should interpret the parameter t to mean the time at 
which the initial incoming field will interact with the system, rather then 
specifying that b{t) is a time-dependent operator at time t. 

0.3.2 Quantum Stochastic Calculus 

The commutator (^l|) is a 5-function because of the Markov approximation. 
The Markovian equations that result have a greatly simplified form, but this 



simplification arises at the expense of fiaving to define stochastic calculus, lead- 
ing to the concepts of Ito and Stratonovich integration much as in the classical 
case. We give here a heuristic review of these ideas. 

Quantum stochastic calculus (QSC) is a non-commutative analogue of Ito's 
stochastic calculus. It was developed originally as a mathematical theory of 
quantum noise in open system and more recently found application to measure- 
ment theory in quantum mechanics. 

For simplicity we consider the situation in which the field is in a vacuum 
state, so that 6(i)|vac) = 0, and thus (6(t)fet(t')) = (vac|6(i)&t(t')|vac) = S{t-t') 
and (&t(i)&(i')) = 0. We define 



B{t):= [ dsbis), B{t)'^ := [ dsb{s)'' 
Jo Jo 



(22) 



For the vacuum averages these definitions lead to 

{[B{t) - B{to)][B{t)^ - i?(to)t]) - |t - tol (23) 

{B{t) - Bito)) - {[B{ty - B(io)t]2) = mt) - Bito)?) = . (24) 
A quantum stochastic calculus of the Ito type, based on the increments 

dB{t) ^ B{t + dt) - B{t), dB{t)'' ^ B{t + dt)^ - B{t)'' (25) 

and dt has been developed by Hudson and Parthasarathy g. The pair B{t), 
B{ty are the non-commutative analogues of complex classical Wiener processes. 

Quantum stochastic integration 

We consider two definitions of quantum stochastic integration: Ito, 



/ f{s)dB{s)= \imy^f{U)[B{U+,)-B{U 



(26) 



and Stratonovich 

(S) / f{.s)dBis)= \iinJ2-[fit,+,) + fitMB{t.+i)~B{tO]. (27) 
Jo "^°°t^2 

In both cases /(t) is a nonanticipating (or adapted) function, i.e. an operator 
valued quantity which depends only on B{s) etc. for s < t. There are analogous 
definitions for the Ito and Stratonovich versions of 

f{s)dB{s)\ f dB{s)f{s),f dBisYfis). (28) 

Jo Jo 

The basic difference between the Ito and Stratonovich form of the integrals is 
that in the Ito form the term f{ti) and [B{ti+i) — B(ti)] arc independent of each 
other, whereas in the Stratonovich form [/(ti+i) + f{ti)] is not independent of 

[Bit,+,) - B{t,)]. 



Example and discussion 

If we use the property d23) , we find 



f dB{s)B{s)^\ =0 



(29) 

(S)^*di?(.)i?(.)t^^i|i|. (30) 

In fact this example shows the three main principles involved. 

i Stratonovich integration follows the rules of conventional calculus. For 
example, the conventional differential of B{t)B{t)^ is the Stratonovich 
differential 

(S) d{B{t)B{t)^) = dB{t)B{t)^ + B{t) dB{t)'' (31) 

so that 



B{t)B{t)^l - (S) / dB{s)B{s)^ + (S) / B{s)dB{s)^ 



(32) 



and using (E3,E^ and taking averages, we find 



\t\-l\t\ + l\t\ (33) 



as should be the case. 



ii Ito Integration increments are independent of and commute with the inte- 
grand. From this independence we see that we can factorize the average 
so that 

£ dBis) BisY^ = I {dB{s)) (i?(.)t) = ^* = (34) 

For Ito differentials the conventional rule (|l]) is replaced by the rules: 

(a) Expand all differentials to second order. 

(b) Use the multiplication rules [for a vacuum state — otherwise use the 
rules (p^ 



dt^ = dB^{t) = dBit)^' = dB\t) dB{t) = 

dB{t)dt = dB{t)'' dt = dt dB{t) = dt dB{t)'' = 

dB{t)dB{t)^ = dt. (35) 



Using these rules, we find 

d[B{t)B^{t)] = dB{t)B{ty +B{t)dB{ty + dB{t)dB{ty 

= dB{t)B{ty +B{t)dB{ty +dt, (36) 

so that using the Ito integral, we get 

(i?(t)S(i)t)lo = |t| = + 0+|i| (37) 



in agreement with (p3). These rules are sufficient for ordinary manip- 
ulation, but are defined only when the state is the vacuum. Rules for 
non-vacuum states will be discussed briefly at the end of this section. 

iii The mean value of an Ito integral is always zero. This follows from the 
independence of dB{t) ~ B{t+dt) — B{t) and f{t), which has been assumed 
to be nonanticipating. 

As in the case of classical stochastic differential equations, an equation 

j^X{t) = F{X{t),t) + G{X{t),t) b{t) + H{X{t), t) b\t) (38) 

cannot be rigorously considered as a differential equation, since the terms in- 
volving bit) and b\t) are in some sense infinite. However this equation can be 
regarded as an integral equation 

X{t) - X{tn) = / dX{t') (39) 

Jta 

{F{X{t'), t') dt' + G{X{t'),t') dB{t') + H{X{t'), t') dB\t')} 

As in the case of non-quantum stochastic differential equations, a simplified 
notation is used, in which the explicit symbols of integration are dropped: 

dX{t) = F{X{t),t)dt + G{X(t),t)dB{t)+H(X(t),t)dB{ty (40) 

which is known as a stochastic differential equation. 

Comparison of Ito and Stratonovich Stochastic Differential Equations 

(SDE) 

The Stratonovich definition is natural for physical situations where the white 
noise approximation is an idealization. While the Stratonovich form has its 
merits of satisfying ordinary calculus, the equation as it stands is rather like 
an implicit algorithm for the solution of a differential equation which makes 
manipulations rather difficult. 

In the Ito definition the differentials dB{t) = B{t + dt) — B{t) etc. point to 
the future. Quantum stochastic Ito calculus obeys simple formal rules which 
can be summarized as follows: 
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i By the CCR (H) nonanticipating functions X{t) commute with the fun- 
damental differentials, [X(t),dB{t)] = etc.. 

ii For two functions X{t) and Y{t) satisfying equations of the type (^), the 
differential of the product X{t)Y{t) is given by 

d[X{t)Y{t)] = [dX{t)] Y{t) + X{t) [dY{t)] + dX{t) dY{t) (4f ) 

where the Ito correction (the last term on the RHS of Eq.(pl|)) has to be 
computed by means of the Ito table 

dB{t) dB{t)^ = dt (vacuum state) (42) 

and all other products involving dB{t), dB{t)'^ and dt vanish. 

0.3.3 Quantum Stochastic Schrodinger, Density Matrix 
and Heisenberg Equations 

The definition of Ito and Stratonovich integration allows us to give meaning to 
a quantum stochastic Schrodinger equation for the time evolution operator U{t) 
of the system interacting with Bose fields 6(i), h\t). The Schrodinger equation 
(EqI) must be interpreted as a Stratonovich SDE 

(S) dU{t) = {-iH dt + dB(i)tc - dB{t)c)}U{t) . (43) 

It is advantageous to convert this Stratonovich equation to Ito form. By con- 
version to Ito form, we can get a modified equation in which the increments 
dB{t), dB{ty and independent of U{t). This Ito form of the quantum stochas- 
tic Schrodinger equation is 

dU{t) = {{-iH-^c^c)dt + dB{t)U~dB{t)c^}U{t) {U{0) = i). (44) 

As the solution U{t) of (H) is nonanticipating, since the increments dB(t), 
dB{ty point to the future, we have [U(t),dB{t)] = [(j{t),dB(t)^ = 0. 

Remark: Conversion from Stratonovich to Ito form: this conversion proceeds anal- 
ogous to the classical case. Consider the Ito equation 

dU(t) = {adt + f3dB{ty - l3UB{ty} U{t) (45) 

We consider an arbitrary Stratonovich integral of a function \(l>{t)) = U{t)\<l){0)) which 
obeys Eq. (^ 

(S) I dB{t)(t>{t) = Urn ^ AS, i [</,(*,+! )+0(t,;)] (46) 

= lim^ AB,[0(iO + ^{aAt,+pABl - (3'' AB,)(j}(t,)] 
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and by AB, ABJ ^ AU 

(S) dB{t)<t,{t) = hm"^^AB,(t>{U) + ^f3AU(t>{U)\ (47) 

dB{t)4>{t) + - / dtP(t){t) 

and similar expressions for the other integrals. Remembering that (ksf) is a shorthand 
for an integral equation we see that (Usl) is equivalent to (ks) if /3 = c, /3^ ~ c^ , 
a — —iH — \c^ c. 

The formal solution of the Eq. (E3) can be written as 

Uit)=Texp[ {-iHds + dB{s)''c~dB{s)c'') (48) 

where T denotes the time-ordered product. This equation gives 
dU{t) = U{t + dt) - U{t) = [exp(-iiJ dt + dB(i)tc - dB{t) c^) - l] lJ{t) 



n=l 



-iHdt + dB{t)'^c-dB{t)c^Yu{t) . (49) 



Applying the Ito rules (^) , terms with n > 2 vanish, the term with n — 1 gives 
-iHdt + dB{t)''c - dB{t)c\ n = 2 gives -^c^'cdt, and we recover Eq. ( ||). 
Obviously, from Eq. ( ^ ) C/(t) is unitary. 

The state vector |*(t)) = f7(i)|*(0)) with initial condition |*(0)) = lip) (g) 
I vac) obeys a stochastic Schrodinger equation 

d|*(t)) = {(-ii? - ic^c) dt + dB{ty c - dB{t) ct}|*(t)) . (50) 

Since |^(0)) contains the vacuum of the electromagnetic field, h{iS)\^{Q)) = 
and thus dB{t)\^{fi)) = 0. Because U{t) commutes with dB(t) it follows that 
dB{t)\^{t)) = 0. This means, as far as dBOt) is concerned, |^(i)) is in the 
vacuum state and therefore the Ito equation (^^ can be simplified to 

d\^{t)) = {-iiJeffdi + dS(i)tc}|*(t)) (51) 

with 

H^s^H-i^c^c. (52) 

This equation has the following physical interpretation. The term dB{ty in- 
volves the incoming radiation field evaluated in the immediate future of t. Thus 
this field is not affected by the system. However, the system does create a self- 
field which causes the radiation damping by reacting back on the system. This 
is the meaning of the term —^c^c\^{t)). The Stratonovich equation does not 
have this term because the evaluation of dB{t)\ half in the future, half in the 
past, itself generates the radiation reaction. 
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Equation of motion for the stochastic density operator 



From the Schrodinger equation (^4|), we can derive the equation of motion for 
a stochastic density operator p{t) = |^(i))(\E'(i)| for the combined system + 
heatbath as 

dp{t) = U(t + dt,t)p(t)JJ{t + dt,t)'' - p(t) (53) 

= -i (H^sp{t) - p{t)Hl^ dt + dB{ty cp{t) cUB{t) 
+dB{tycp{t) + p{t) dB{t) c^ . 

Here 

U{t,to) = U{t)U{t^)^ (54) 



is the time evolution operator from ip to t. Eq. (53) is derived by expanding 
U{t,t+dt) to second order in the increments (compare Eq. (Q)), and using the 
Ito rules (|4^). If we now trace over the bath variables, this will perform an 
average over the dB{t), dB{ty operators. We use the cyclic properties of the 
trace, and derive the usual master equation for p{t) = Trsipit)}, 

j^p{t) = -i(H,iip-pHl^^+Jp (55) 

= Cp 
with C a Liouville and "recycling operator" 

Jp = cpc'^ . (56) 

For a system operator a we define a Heisenberg operator a{t) = tj{ty a U{t) 
that obeys the Ito quantum Langevin equation 

da[t) = (j{t+dt,t)'^ a{t)U{t+dt,t)-a{t) (57) 

= -i[a,Hdt + idB{t)c^ - idB{t)^c\ + {c^ac~ a-c^c c^ca)dt . 

Remarks 

(i) Generalization to more channels: For reference in later sections we need the 
master equation for Nc output channels. It is given by 



pit) = -i[H,p\+Y.^c,pc]~-c]cjp~p-c]cj) (58) 

= -{H,sp-pHl^)+Y.^iP = ^P (59) 



j=i 
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with 'effective' (non-hermitian) Hamiltonian 

1 "^ 
H,s^H-i-Y,c]c,, (60) 

and 'recycling operators' 

JoP^Cjpc] (j = l,...,7Ve). (61) 

A master equation of this form is derived by couphng a system to Nc reser- 
voirs with system operators Cj {j — 1, . . . ,Nc)- 
(ii) Non-vacuum initial states: So far we have assumed a vacuum state (a pure 
state) in the input field. The above formalism can be generalized to nonpure 
initial states, provided we stay in a white noise situation |Q, |2| Q. For a 
Gaussian field with zero mean values {b{t)) = (&(t)^) = 0, and correlation 
functions {b{t)^b{f)) = N6{t - t'), {h{t)b{t')'^) = {N + l)5{t ~ t'), and 
{b{t)b{t')) = MS{t — t') where N is the mean photon number and M is a 
squeezing parameter (|Mp < N{N + 1)), it can be shown [0, |l^, ^ that 
the Ito table has to be extended according to 

dB(t)UB(t) = Ndt dB{t)dB{t)'^ = iN+ l)dt 

dB{t) dB{t) = Mdt dB{t)'^ dB{ty = M*dt (62) 

dB{t) dt^dB^t) =0 . 

0.3.4 Number Processes and Photon Counting 

The description of photon counting requires an extension of quantum stochastic 
calculus to bring in the so-called gauge or number processes. The number 
process arises by noting that the total number of photons counted between time 
and t is 

A{t) := [ dsb''{s)b{s) . (63) 

Jo 

This assumes a perfect photodetector with unit efficiency and instantaneous i5- 
function response. The operator A{t) leads to the construction of the quantum 
analogues of a Poisson process. Physically, dA{t) is an operator whose eigen- 
values are the number of photons counted in the time interval {t,t + dt]. It is 
only possible to make sense of dA{t) in situations in which the initial state is a 
vacuum. The reason for this is physically quite understandable as a white noise 
field would lead to an infinite number of counts in any finite time interval. 
We will also need Ito rules for dA{t) [lO[ |, 

dA{t)dA{t) = dA{t), dB{t)dA{t)^dB{t) (64) 

dA{t)dB[t)^ = dB{t)'' (65) 
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together with (E2|) and all the other products involving dB{t), dB(ty , dA{t) 
and dt vanish. 

We can deduce these Ito rules by noting that for any optical field with no white 
noise component, the probability of finding more than one photon in the time interval 
dt will go to zero at least as fast as dt^. We can use this fact to show that if we 
discretize the integral \A) — J[dA{t)'^ — dA{t)]\(l>{t)) , then in the limit of infinitely 
fine discretization the mean of the norm of |^) goes to zero. This means that in 
a mean-norm topology, \A) — > 0, and thus we can formally write inside integrals 
dA(t)^ — dA{t). This equation ensures that the only eigenvalues of dA{t) are and 1, 
and that we can only count either one or no photons in a time interval dt. 



0.3.5 Input and Output 



In quantum optics the Bose fields (20) plays the role of the electromagnetic 



field. These fields represents the input and output channels through which the 
system interacts with the outside world. In particular, b{t) = b™{t) is the input 
field which is the field before the interaction with the system. In the work of 
Gardiner and CoUett W O , and Barchielli M, M, nrt the Heisenberg equations 
for system operators are expressed as operator Langevin equations with driving 
terms b{t). Gardiner and CoUett also define output operators for the fields after 
the interaction with the system, 

1 f^ 
fo™t(t) :^ / dw6(w,T)e-^('^-'^°)(*-'^) (66) 

V 27r J-oo 

with b(uj, T) the Heisenberg operator of 6(w) evaluated at some T in the distant 
future. In Ref. |l^ it is further shown that 

6°"t(i)=6'"(t)+c(t) (67) 

which expresses the out field as the sum of the in field and the field radiated by 
the source. The fundamental problem of quantum optics consists in calculating 
the statistics of the out field for a system of interest, in particular the normally 
order field correlation functions 

(6°"*(ii)t . . . 6°"*(i„)t6-t(i„+i) . . . 6°"*(i„+™)) , (68) 

given (all) correlation functions for the in field. 
Formally, we introduce output processes [ |lO| 

B°"*(t) := C/(t)t B{t) U{t), A°"'(i) := U{t)^ A{t) U{t), (69) 

and the output fields by formal (forward) derivatives 

6°"t(i) := dB°''\t)/dt = lim fB°"Hi + h) - B°''Ht)]/h . (70) 
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Similar definitions hold for B°^*(i)t and &°^*(i)^ From Eqs. (|9|,|4|) we see tliat 
tfie time evolution operator U{t, s) from time s to time t only depends depends 
only on the fields b{r), 6(r)' with r between s and t. By the commutation 
relations 

[U{t, s),B{t)] =0 for r < s (71) 

this gives 

B°"t(t) = U{t)^U{T,t)^B{t)U{T,t)Uit) 

= U{T)^B{t)tj{T) {T>t). (72) 

A QSDE for the output processes is derived using QSC \\S\ : 

dB°''\t) = B°''\t + dt) - B'^'^t) (73) 

= U{t+dt)'' B{t + dt)U{t+dt) - U{tyB{t)U{t) 

= U{t)''U{t+dt, t)UB{t)U{t+dt, t)U{t) 

= dB{t)+c{t)dt (74) 

where the second line follows from (69), the third line from (|7l|), and the last 
line is obtained from (|l^) and the Ito rules (^2|) in agreement with Eq. (|6j). In 
a similar way one finds 

(75) 
dA™*(t) = dK{t) + c{t) dB{t)^ + c(t)t dB{t) + c{tyc{t) dt . (76) 

which expresses the operator for the count rate for the out field as the sum of 
the in operator, an interference term between the source and in field, and the 
direct system term. 

0.4 Counting and Diffusion Processes: Stochas- 
tic Schrodinger Equation and Wave Func- 
tion Simulation 

0.4.1 Photon Counting and Exclusive Probability Densi- 
ties 

In this section we discuss photon counting of the output field as realized by a 
photoelectron counter ^ |[ (see Fig. ga). Mathematically, this corresponds to 
a class of continuous measurements known as "counting processes" [ri|. We 
consider a system which is coupled to j = 1, . . . , Nc output channels with coun- 
ters which act continuously to register arrival of photons. A complete descrip- 
tion of counting process is given by the family of exclusive probability densities 
(EPDs): 
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Figure 2: (a) Photon counting, and (b) homodyne detection. 



P(* (0|p) = probability of having no count in the time interval {to, t] when 
the system is prepared in the state p at time to- 

Pto iJrn 7 im,; • • ■ ; ji , ii |p) — multitime probability density of having a count 
in detector ji at time ti, a count in detector J2 at time i2 etc. and no 
other counts in the rest of the time interval (to, t] (with jk = 1, ■ ■ • , Nc] 
to<ti < ... <t^<t). 

Knowledge of EPDs allows the reconstruction of the whole counting statistics. 

Mandel's counting formula 

Physically, the operator for the output current of a photodetector is proportional 
to the photon flux. To simplify notation we restrict ourselves in the following 
to the case of a single output channel Nc = 1. Furthermore, since the aim of 
the present discussion is to discuss the relation with the stochastic Schrodinger 
equation, we will concentrate on the case of ideal photodetection with unit 
efficiency. The photon flux operator is 



dA°"*(t) = 6°"*(t)t6°"'(t) dt =: /(t) dt 



(77) 



Here A°"*(t) is the operator corresponding to the number of photons up to time 
t in the out field. We note that the operators A°"*(t) are a family of commuting 
selfadjoint operators. Photoelectric detection realizes a measurement of the 
compatible observables A°"'(t). In quantum optics the usual starting point is 
the photon counting formulas first derived by Kelley and Kleiner, Glauber and 
Mandel 0: 



Pmit, to) — 



A( /^ dsl{s) 



exp 



ds I{s] 



to 



(78) 
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gives the probability for m photoncounts in the time interval (io , t) where : . . . : 
indicates normal ordering, and (...) = Trs^sj. . . po} with pQ — p® |vac) (vac| is 
the quantum expectation value with respect to the state of the system. Eq. ( [Tq ) 
is valid for a detector with unit efhciency. In particular, the probability of no 
count in a time interval (io — 0, t] is 

^t*o=o(0|p) = /: cxp (- f dsiis)) :\ , (79) 







and the multitime probability density for detecting (exactly) m photons at times 
ti, . . . , tm in the time interval (to — 0, t] is 

pl=o(tm,---,ti\p) = (-.i{h)...i{t„,)exp(-l dsi{s)) :\ . (80) 

This establishes the relation to the EPDs introduced above. 

On the other hand, the non-exclusive multitime probability density for de- 
tecting a first photon at time ii, a second (not necessarily the next) photon at 
time ^2 stc. and the m-th photon at t^ is given by the intensity correlation 
function 

' i{t,)...i{t„,):) . (81) 



The characteristic functional and system averages 

As a next step we wish to express the out field correlation functions in terms 
of system averages. The standard approach is to express in Eqs. (u9t) and ( pO[ ) 
the out field in terms of the in field plus the source contribution, 6°"'(i) = 
6'"(t) + c(t), and to prove and apply the quantum fluctuation regression theo- 
rem (QFRT) to evaluate the normally- and time-ordered multitime correlation 
function for the source operator c(t)j£|. Instead we will follow here the elegant 
procedure outlined by Barchielli [y, |lC | of defining and evaluating a character- 
istic functional for the EPDs. 

We define a characteristic functional as the expectation value of the charac- 
teristic operator V°^^ [k] according to 

$t[fc] = (yr^[k]\ with Vr*'{k]:=exp(i f k{s)I{s)dsj . (82) 

When we use the normal ordering relation 0] 

exp (i I k{s) i{s) dsj =: exp ( / (e''^^") - 1) /(s) dsj : , (83) 

expand the exponential, and use expressions (|79|) and (80) we obtain 



$t[fc] :=Fo*(0|/o)+^ /dim/ dtm-i ■ ■ . dti 
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expii^k{tr)\pl{t„„...,ti\p) (84) 

which relates the characteristic functional to the exclusive probability densities. 
Thus ^t[k] determines uniquely the whole counting process up to time t. Func- 
tional differentiation of $t[fc] with respect to the test functions k{t) gives the 
moments, correlation functions and counting distributions. 

According to Eq. (pOh the out and in processes are related by A°"*(i) = 
U (t)^ A{t)U (t) and the characteristic operator satisfies the QSDE 

dVr'[k] = Vr'M U'"^*'^ - l) dA°''\t) (Fo°"*[/c] = i), (85) 

dA°^*(t) = dA{t) + c{t) dB{t)^ + c(i)t dB{t) + c{t)'^c{t) dt . (86) 

In a similar way we can express the characteristic operator for the out field in 
terms of the in field; 

Vr'[k] = U{t)Uxp(i f k{s)dA{s)]u{t) (87) 

= U{t)^Vt[k]U{t) . 

This allows us to write 

Mk] = TrseB{Vt[k]U{t)p{Q)U{t)^ (88) 

= Trsxtik] with xt[k]:=Ti-B{Vt[k]p{t)}. (89) 

which expresses the characteristic functional as the system trace over a char- 
acteristic density-like operator Xt[k] in the system space. Using QSC we can 
derive an equation for Vt[k] p{t), 

d[Vt [k] m] = dVt [k] pit) + Vt [k] dp{t) + dVt [k] dp{t) (90) 

with dp{t) from (|53|), and Vt[k] obeys an equation similar to (p5|). Taking the 
trace over the bath, Tr^j. . .}, as in the derivation of the reduced density matrix 
(p5|), we obtain the following equation for the characteristic reduced density 
operator 

j^Xt[k] = [C + (e"=(*) - 1)^) xt[k] (xo[fc] = P) , (91) 

with L the Liouville operator from Eq. (|55|), and J the "recycling operator" 
(p6[). Note that for k(t) =0 this characteristic density operator coincides with 
the system density operator, xt[k = 0] = p(f), and Eq. ( |9l| ) reduces to the 
master equation (pq). 

We solve Eq. J9l|) by iteration 



Xt[fc] =5(t,0)p-|- V] / rfi,„ / dt^^i... diiexp iV'fc(t 

m=l -^0 -^0 -^0 V r=l 



Sit, tm)J{tm)S{trn, t,„-l) . . . Siti, ti)J {ti)S{ti, 0)p 
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with propagators 

S{t, h)p = U,s{t, to)pUlg{t, io) (i > to), (92) 

Ucs{t,to) ^ -iHesit)Uesit,to) (C/off(io,to) = i) . (93) 

Taking a trace over the system degrees of freedom, and comparing with (p3) we 
can express the EPDs in terms of system averages 

P^{0\p) = Ti-s{S{t,0)p}, (94) 

pl{ti,t2,...,tm\p) = 

TTs{S{t, t^)J{t^)Sitrn,t^-l) . . . S{h,h)J{ti)Sih,0)p} . 

The structure of these expressions agrees with what is expected from continuous 
measurement theory according to Srinivas and Davies pq ] . 

Generalization to many channels 

As a reference for the foUowing sections, we give the generaUzation of these 
equations to Nc channels: 

P4^o(0|p) = Trs{5(t,0)p}, (95) 

p*o=o(ji'^i;---;im,<m|p) = 

Tr5 {Sit, t^)Jj^ (t,^)S{t^,t,^^i) . . . S{t2M)Jj, (ii)5(ti, 0)p} 
where io = < ti < . . . t„i < t and j^ = 1, . . . , Nc- 

Examples of the use of the characteristic functional 

We conclude our discussion with a few examples for the application of the char- 
acteristic functional: The mean intensity of the outgoing field is 

^(*^) = -'-S^^'^T^''^\^=o^^'s{Jp{t)} (T^t). (96) 

The intensity correlation function can be written as the sum of a normally 
ordered contribution plus a shot noise term 



mi{t')) = (:/(i)/(t'):)+<5(t-t')(/(i)) 



= Tr{Jit)T(t, t')J{t')p{t')} + 5{t - t')Tvs{J{t)p{t)} . 
The normally ordered m-th order intensity correlation function is 

(: i{ti) . . . i{tra) :) = Trs{ J(t„,)T(t„„ <,„_i) . . . T{t2M)J(ti)T{ti,Q)p} 

(98) 
with T(i, to) = exp>C(< — t^) the time evolution operator for the density matrix 
according to the master equation. 
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0.4.2 Conditional Dynamics and A Posteriori States 

If during a continuous measurement a certain trajectory of a measured observ- 
able is registered, the state of the system conditional to this information is 
called a conditional state or a posteriori state denoted by pdt), and the corre- 
sponding dynamics is referred to as conditional dynamics [^, llj 1^, |lj]. In the 



Pio,{t,t+i]\p) ^ "^"t.ry ;. ;Tr (99) 



present case of photodetection (counting processes) we assume that the coun- 
ters have registered in the time interval (to = 0, t] the sequence ji, ii; . . . ; j„, t„ 
{ti < . . . < tn) ■ To derive an expression for Pc{t) Barchielli and Belavkin pi[ 
consider 

pI)^* {ji,ti;---;jn,tn\p) 

Pl{il,h\---]jn,tn\p) 

which is the conditional probability of observing no count in {t, t + t] given p at 
time to = and having observed the sequence ji, ti; . . . ; j„, in- Obviously one 
can write P(0, (t, t + tJlp) = Pt*+*(0|pc(t)) with 

Pc{t) = S{t, tn)Jj„ (t„)5(t„, t„_i) . . . S{t2,ti)Jj, {ti)S{ti,0)p/Tis{- ■ ■} (100) 

where Trgj. . .} is a normalization factor to ensure TrsPc(t) == 1- Similar argu- 
ments can be given for other EPDs. Thus we identify ( |100| ) with the a posteriori 
state. 

Discussion g |l^, |l|, [l4| : 

(i) The interpretation of Eq. (|lOO| ) is as follows: When no count is registered 
the system evolution is given by 5(t, to) and the state of the system between 
two counts is 

Pc{t) = S{t, tr)pc{tr)/Trs{S{t, tr)pc{tr)} (t > t,) (101) 

where t^ is the time of the last count, and Pc{tr) is the state of the system 
just after this count. 

When a count j at time t = t^ is registered the action on the system is given 
by the operator Jj , and the state of the system immediately after is 

p,{tr + dt) = Jj{tr)pc{tr)/^rs{Jj(tr)Pc{tr)} ■ (102) 

We will call this a "quantum jump." 
(ii) When the system is initially in a pure state described by the state vector f/;, 
p — |V')(^|, it will remain in a pure state ipc{t)i Pc{t) = |^c(0)(V-'c(i)|- The 
time evolution between the counts is 

Mt) = u,ff{t, tr)Mtr), Mt) = Mt)/\\Mt)\\ (103) 

i.e. 

^Jcit) = -iH,sMt) (104) 

dt 
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where tpcit) and ipc{t) denotes the unnormaUzed and normaUzed wave func- 
tions, respectively, and a count j, U is associated with the quantum jump 

lpc{tr + dt) = XjCjXpdtr) ■ (105) 

The complex numbers Xj 7^ (with dimension [time]^'^) are arbitray and 
can be chosen, for example, to renormalize the wave function tpcit) after the 
jump, 
(iii) In preparation for the wavefunction simulation to be discussed in Sec. |0.4.3 



we note: The mean number of counts of type j in the time interval (t, t+dt] 
conditional upon the observed trajectory ji, ti; . . . ; j„, t„ is 

.!«^...i.,(o)*={^;««';<''>''' . (106) 

and the joint probability that no jump occured in the time interval (0, i], 
and a jump of type j occured in (i, t + dt] given p at time t = is 

pro. .1 ,s) ) X PM„ = { ^;£ «5';,ii?i;' * ■ im 

with pc{t) ^ S{t,0)p/Trs{S{t,0)p}. 

0.4.3 Stochastic Schrodinger Equation and Wave Function 
Simulation for Counting Processes 

Stochastic Schrodinger Equation 

We can combine the time evolution for the (unnormalized) conditional wave 



function '(/'c(i) according to Eqs. (103,104) into a single (c-number) Stochastic 



Schrodinger Equation (SSE) jTll, |14|. A typical trajectory Nj{t), where Nj{t) is 
the number of counts up to time i, j = 1, . . . , N^ is a sequence of step functions, 
such that Nj (t) increases by one if there is a count of type j and Nj (t) is constant 
otherwise. Therefore, the Ito differential 

dN.it) . N,it + dt) - N,it) = { I l^^^ ^ - (^'^+^^1 (108) 

fulfills [dNj{t)]'^ — dNj{t). Moreover the probability of more than one count in 
an interval dt vanishes faster than dt. We thus have the Ito table 

dNjit) dNkit) = 6jk dNj{t), dNj(t) dt = (109) 

and obtain a Stochastic Schrodinger Equation for the unnormalized wave func- 
tion Ipcit) 

dMt) = Mt+dt)~Mt) (110) 

= I -iH^a dt + J2 (-^jCj- - i) dNj (t) i 4>c{t) 
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with Xj 7^ 0. According to (106) the mean number of counts of type j in the 
interval {t,t+dt] conditional to the trajectory ji,ii; .. .■,jn,tn (indicated by the 
subscript c) is 



{dNj{t))=\\cj^Pc 



dt . 



(Ill) 



The SSE ( |110[ ) can be converted to an equation for the normalized ipdt) = 
'!/'c(i)/||V'c(i)||- As an illustration for Ito calculus we give details in Appendix [l|. 
The result is 



dMt) 



^c{t + dt) - Vc(i) 



N, 



-iHdt- 






^cj- 



4^^- 



""h 



dt 



(112) 



1 dN. 



iW}V'c(i) 



where 



= {Mt)\---\Mt)) 



(113) 



Eq. (112) is a nonlinear Schrodinger equation. (Note: Eq. (112) is not unique 
up to a phase factor: ipc(t) -^ e*'^(*)-0c(i) where (f>{t) is solution of a stochastic 
differential equation d(j){t) — a{t) dt + ^bj{t)dNj{t) with arbitrary a{t) and 
b,(t).) 



Equation of motion for the Stochastic Density Matrix 

Finally, we derive an equation for the Stochastic Density Matrix pc(t) — |'0c(i))(V'c(O|j 



dpcit) = pc{t+dt) - pc{t) 
= Cpcit)dt + J2 



(114) 



JjPcit) 



fr{ VTrsUPc(t)} 



-pdt)]{dN,(ty 



TTs{JjPc{t)}dt) . 

Discussion: The Master Equation (|5q ) is (re)derived by taking the stochastic 
mean of Eq. {114) ■ We first take the mean for the time step t -^ t+dt conditional 



upon a given trajectory. All quantities on the right hand side of Eq. (114) 
depend only on the past (they are nonanticipating or adapted functions), but 
the a posteriori mean value of {dNj{t)) is given by (111), and thus the last term 
vanishes. Then we take mean values also on the past, and obtain the master 
equation (pq): p{t) = Cp{t). Thus, if the results of a measurement are not read, 
i.e. no selection is made, the state of the system at time t will be p{t) — {pc{t))g^, 
and p{t) is the a priori state for the case of continuous measurement. 

From our construction it is obvious that the statistics of the jumps Nj{t) is 
identical to the statistics of photon counts discussed in Sec. p. 4.1. However, as 
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a consistency check and for pedagogical purposes, we will rederive the count- 
ing statistics in the framework of the characteristic functional for the counting 
process Nj{t) and show that it is identical to the characteristic functional de- 
rived in Sec. 0.4.1. We confine ourselves to Nc = 1. In the present context the 
charcteristic functional for the counting process Ni{t) = N{t) is defined as 



<^t[k]^{Vt[k])^,with Vt[k]:^exp 



I I k{s)dN{s) 





(115) 



which should be compared with the definition ( p2| ) for the output operator 
A°"* (i) . The quantity Vt [k] satisfies the c-number SDE 



dVt[k] = Vt[k] U^^'^ - l) dN{t) . 
In analogy to Eq. (Isl) we introduce 

$t[fc]-Trsx4fc] with xt[k]:^{Vt[k]p,it))^^ 
According to Ito calculus 



(116) 



(117) 



d[Vt[k] p,{t)] - dVt[k] p,{t) + Vt[k] dp,{t)] + dVt[k] dpc{t) (118) 



where (iVt[fc] and dpc{t) obey Eqs. (116) and (114), respectively. Repeating 
the arguments given in the derivation of the density matrix equation from 
the stochastic density matrix equation given above, we take a stochastic av- 
erage (■ . . )„t of the equation for Vt[k] pc{t) and find that Xt[k] = iyt\k\ Pc(0)st 
(Eq. (117)) satisfies the equation of motion and initial condition (|9l| ) given in 
Sec. 3.4.1 . This completes the proof that the characteristic functional defined 
in Eq. (|ll5| ) is identical to the characteristic functional of Sec. |0.4.1 . 



A direct derivation of the c-number SSE (112) from the QSSE (51) has been 
given by Gardiner et al. [|l3| . Goetsch and Graham iQ have given a derivation 
by working with a eigenrepresentation of the set of commuting operators dA(t). 



Wave Function Simulation 



The conditional dynamics defined by Eqs. (103,105) together with ( |107| ) suggests 
a wave function simulation of the reduced system density matrix p{t) as follows 
§ 0, |n|, [||14[ |l|, [8Kl9[ §, H, H g, H, H, H, §, & |lH|, H : 

Step A: We choose a system wave function ipdto), initially normalized ||?/'c(io)|| = 
1, and set the counter n for the number of quantum jumps equal to zero, 
n = 0. 

Step B: We propagate the (unnormalized) conditional wavefunction ^c accord- 
ing to Eq. ( |103| ) and simulate the time t and type of count j of the next 
quantum jump according to the conditional density (107), 



cU,t) = Wcj'ipc 



{t>tn,\\Mtn)\\ = l)- 



(119) 
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Note: One possible way to determine t and j is to proceed in two steps. First, 
we find a decay time t according to the delay function c{t) = Yl^iJ'^)- This 
is conveniently done by drawing a random number < r < 1 from a uniform 
distribution and monitoring the norm of V'c(i) until 



t 

'~/j.l\ _ 1 II 7 ^j.M|2 



dtc{t') = l-\\il)c{t)f (120) 

= '■<:[0,1] (t>in; ||^c(in)i| = 1). 

Second, the type of count j is determined from the conditional density c{j\t) = 
c(j, t)/^j. c{k,t) for the given time t. 

Incrementing n — > n + 1 we identify i„ = t and j„ = j with the decay 
time and type of cou nt, re spectively. The wave function after the quantum 



jump is give n by Eq. ( 105 ). We renormahze this wavefunction and continue 
integrating ( 103 ) up to the next jump time, i.e. return to the beginning of 
Step B. 
Step C: An approximation for the system density matrix is obtained by repeat- 
ing these simulations in steps A and B to obtain 

Pit) = {\Mt)){Mmst (121) 

where (. . .)^^ denotes an average over the different realizations of system 
wave functions. 

Remarks: 

(i) A simulation of the quantum master equation in terms of system wave func- 
tions can replace the solution of the master equation for the density matrix, 
and an important feature of the wave function approach, first emphasized 
by Dalibard, Castin and M0lmer fl^, is that one has only to deal with 
a wave function of dimension N, as opposed to working with the density 
matrix and its 7V^ elements. Thus, simulations can provide solutions when 
a direct solution of the master equation is impractical because of the large 
dimension of the system space. Convergence of the simulation vs. a direct 
solution of the master equation in terms of "global" and "local" observables 
are given in Ref. |18| . 

In many problems one is interested in a steady state density matrix. In 
this case it is often convenient to replace the ensemble averages by the time 
average of a single trajectory. 

A simulation has the further advantage that it allows the calculations to 
be performed on a distributed system of networked computers ||25| , with 
a corresponding gain in computational power. In view of the statistical 
independence of the different wavefunction realizations, parallelization of 
the algorithm is trivial. 
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(ii) The iVc-channel master equation l5q) is form invariant under the transfor- 
mation 

Cj^J2^jkCk where (J^UjkUtk = 6,i) (122) 

k=i j=i 



(that is with U a unitary matrix). The decomposition of Eq. (pq) to form 
quantum trajectories ipdt) is thus not unique. Different sets of jumps oper- 
ators {cj} not only lead to a different physical interpretation of trajectories, 
but an appropriate choice of Cj may be crucial for the formulation of an 
efficient simulation method for estimating the ensemble distribution. An 



example will be given in Sec. 0.5.2 |W3] (see also |18|) 



(iii) Any master equation conforming to the requirements for p is of the Lind- 
blad form H] 



N, 



[H, P]+Y. ^-1 (^C'/'c] - c]c,p - pc]c^ . (123) 



This can be brought into the form ( p8| ) by diagonalizing the Hermitian 
matrix 7^ by a unitary transformation V^^, 

7»j=5I^^7'«7C^70^ (i,J- =l,...,iV,) (124) 

7=1 

with eigenvalues K-y > 0, and defining 

i=l 

Applications of this procedure to the master equation for squeezed noise can 
be found in Ref. [|o). 

0.4.4 Simulation of Stationary Two-Time Correlation Func- 
tions 

We are often interested in system corr elation functions {A{t)B{tQ)) and their 



Fourier transform (spectra) (see section 3.5.2 ). According to the quantum flue 



tuation regression theorem the system correlation function can be written as 

(A(t)B(io)) = Tr5{Ae^(*-*«)i?p(io)} {t > io) • (126) 

This correlation function can be obtained by solving the density matrix equation 
for pit) and an equation for the first order response p'+^(i) 

m = Cp{t), (127) 

p^+\t)=Cp^+\r) + BJ{t)p{t). (128) 
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By formally integrating Eq. (128), 

p(+)(t)- / dhe^^'-'^^Bf{h)pih) 



and choosing a ^-function f{t) = 5{t — to) we have 

Trs{V+nO} = (^WS(to)) it>to) 



(129) 



(130) 



We will show that these correlation functions can be obtained by solving 
(simulating) the set of stochastic Schrodinger equations O, p^, |2^ 



dMt) = 



-iHcs dt + ^{XjCj - i)dNj{t) 



j=i 



iHcS dt + ^{XjCj ~ i)dNj{t) 



Mt) (131) 

4>i+\t) + f{t)Bi;,{t)dt 



with quantum jumps dictated by 4'c{t) according to {dNj{t)) — \\cj'ipc{t)\\'^ dt. 
For a ^-kick tpi (t) obeys the same Schrodinger equation as V'c(i) (Eq. (131)), 



where the inhomogeneous term in (132) translates into the initial condition 

4>i+\to + dt) = BMto) ■ (132) 

It is straightforward to show that the stochastic averages 

p{t)^(\Mt))mt)\/\\Mt)V)^^ (133) 

p^+Ht) = (Hi+Ht)){Mt)\/\mt)r) (134) 



(both normalized with respect to ipc{t)) obey Eqs. (127, |128| ), and thus 

{Ait)Bito)) = { {{Mt)mi-'Ht))/\\Mt)r) ),„ (135) 



where 



indicates an average over radmonly selected initial "kick" -times 



to- A physical interpretation of this procedure is summarized in the discussion 
following Eq. ( |170| ). The above derivation emphasizes the simulation of a corre- 
lation function. We can simulate a spectrum directly when instead of the J-kick 
we use a function f{t) ex exp {—ivt) [|l3[ |2C|, pl| , [25[ |. For another approach to 
simulate correlation functions we refer to |18[ . 

0.4.5 Diffusion Processes and Homodyne Detection 

Le us consider homodyne detection as shown in Fig. Ep P, nO, n4|. In the 
simplest configuration, the output field of the system 6°"* (i) is sent through a 
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beam splitter of transmittance close to one. The other port of the beam splitter 
is a strong coherent field which serves as a local oscillator. For a single output 
channel, the transmitted field is then represented by 6°"*(i)' = 6°"'(i) + 7 = 
c + 7 + b{t) . Assuming a real field 7 the operator for photon counting of the 
outgoing field is 

dK°''\t)'/dt = IP^'^tY IP^'^t)' (136) 

= 7' + 7 [c + c^ + h{t) + fe(t)t] + [ct + 6(i)t] [c + h{t)\ . 

The ideal limit of homodyne detection is for infinite amplitude of the local 
oscillator. Physically, in this limit the count rate of the photodetectors will go 
to infinity, but we can be define an operator for the homodyne current /''°™(i) 
by subtracting the local oscillator contribution: 

r-°'-\t) = hm "^ ^^ , ' ^ . (137) 

7^00 7 dt 

Defining quadrature operators for the in and out field, 

dX{t) := dB{t) + dBt(t), dY{t) := ~idB{t) + idB{t)'^ (138) 

and for the system dipole 

x:=c + c\ y:=—ic + ic*, (139) 

we can write 

dl°"'(t) := I^°''^{t)dt (140) 

= x{t)dt + dX{t) 

where the second line is the QSDE for the X-quadrature of the out-field which 
follows from Eq. ([7^. This equation shows that homodyne detection with 7 real 
corresponds to a measurement of the compatible observables 6°^*(i) + 5°"*(i)t = 
dX°^^{t)/dt (t > 0), the quadrature components of the field. 

The complete statistics of the homodyne current is obtained from a charac- 
teristic functional pOl 

$t[/c] = (vr'[k]) with Vr'[k]-=exp(i f k{s)dX°'''{s)) .{Ul) 

The characteristic operator is determined by the QSDE ^ 

dVr'[k] = Vr'[k] (ik{t) dX°'-'\t) - ifc(i)2 dt] (T>o™*[fc] = i)(142) 



^ We note the connection between Eq. ( 142 ) and normal ordering relation 



exp [ i / fc(s)dX°"*(s) 1 =: exp ( i / fc(s)dX°"'(s) ~- k(sf ds 
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together with Eq. (14C) for dX°"*(i). FoUowing the derivations of Sec. D.4.1 



we define a characteristic density operator Xt[k] which gives the characteristic 



functional $f[fc] (141) according to (p8[). We obtain the equation 



-Xt[k] = Cxt[k] - ^k{t)\t[k] + ik{t) {cxt[k] + Xt[k]c^) where {xo[k] = p(0)) 

(143) 
corresponding to a Gaussian diffusive measurement |ll| . An expression for the 
mean homodyne current and current - current correlation function will be given 



below in Eqs. (159) and (|160D, respectively. 



0.4.6 Stochastic Schrodinger Equation and Wave Function 
Simulation for Diffusion Processes 

We now turn to a c-number stochastic description and the derivation of a c- 
number stochastic Schrodinger equation for diffusion processes. Such an equa- 
tion was first derived by Carmichael [gj (see also ||lj|) from an analysis of homo- 
dyne detection, and independently in a more formal context by Barchielli and 
Belavkin p|. 

From ( |106D the mean rate of photon counts is 

(dNit))^ = Trs{ic^+j){c + j)p,it)}dt (144) 

= Trs{{j^+-fx + c^c)pcit)}dt, 

and registration of a count in homodyne detection is associated with a quantum 
jump V'c ■*— [c + 7)V'c of the system wave function. Furthermore, we note that 
the master equation is invariant under the transformation 

c -^ c + 7 (145) 

H -^ H -i-{j*c--fc') (146) 

which involves a displacement of the jump operators by a complex number 7. 



With these replacments we obtain from the SSE (112) 



dMt) = (-iH-^{c^c + 2^c-{c^c)^^j{x)j) dti;,{t) (147) 



""^^ l]dN{t)Mt) 



We are again interested in the ideal limit when the local oscillator amplitude 
7 goes to infinity. As pointed out before, in the limit that 7 is much larger than 



c, the count rate (144) consists of a large constant term (the counting rate form 



the local oscillator), a term proportional to the x quadrature of the system 
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dipole, and a small term, the direct count rate from the system. Subtracting 
from the counting rate the contribution of the local oscillator we write 

dN{t)=-f^dt + -fdX{t) (148) 

which defines the stochastic process X{t). For 7 -^ cxd the stochastic process 
X{t) has the Gaussian properties 

(i) {dX{t))^ - {x)^ dt (149) 

(ii) [dX{t)f ^ dt, dX{t)dt^O (150) 

=> dX{t) - {dX{t))^ = dW{t) (151) 

where dW{t) is a Wiener increment [dW{t)f = dt, {dW{t))^ = 0. Thus the 



jumps (147) are replaced by a diffusive evolution. Th e deri vatio n of (149) is 
straightforward [|ll: Eq. ( |l49| ) follows from Eq. ( |l44| ); Eq. ( |l50|) is obtained 
from the Ito table ([l09| ), [dX{t)Y = dN{t)/-f^ = dt + dX{t)/j. Physically we 
interpret 

/-(O = ^^ . hm !m:Ulj^ (152) 

^ ^ dt 7^00 ■J dt ^ ^ 

as a stochastic homodyne current with ^(t) = dW{t)/dt white noise. 



Taking this limit in Eq. (147), we obtain a SSE for diffusive processes 



dMt) = {(^-iH-^ic^c-{x)^c + ^{x)l)^ dt (153) 

+ ic-{x/2)^)dWit)}ij,{t). 

Eq. ( |153| ) is not unique; in particular we can make a phase transformation 
V'c(i) -^ e*'^^*Vc(0 which allows us to rewrite this equation as IT^| 

dMt) = {(^-iH-^c^c + {c^)^c + ^\{c)J'^ dt (154) 

+ {c^{c)^)dW{t)}Mt), 



and there are also versions with complex noises l34, |35, 37 . 

The unnormalized version of the equation has been given by Carmichael B , 
and more recently by Wiseman ||lj, |l^ 

di>c{t) = -iH - -c^c + /''^(i) c] dt Mt) (155) 

which demonstrates clearly how the state is conditioned on the measured pho- 
tocurrent I^°'^{t). 
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A stochastic density matrix equation for diffusive processes is directly ob- 
tained from these equations, 



dpcit) = Cpc dt + [(C - (c) J Pc + Pc (ct - (ct) J] dW{t) 



(156) 



and taking stochastic averages we see - as expected - that the a priori dynamics 
satisfies the master equation (|55|). 

As shown by Goetsch et al. jl^l the above c-number SSE for diffusive 
processes can be derived directly from the QSSE (|5l|). The idea is to replace 
dB{t)U\^{t)) -^ (dB(t) +dB(t)t)c|^(i)) = dl(i)T7|*(i)) in Eq. (|5l|) and to 
subsequently project on an cigcnbasis of the set of commuting operators dX. 
Generalization of the SSEs to a squeezed bath have been given in ||lj, Q. 

The stochastic homodyne current /''°™(i) ( |153| ) has the same statistical prop- 
erties as the operator version I^°"^{t) as becomes evident by defining a charac- 
teristic functional fOl 



I I k{s)dX{s) 



*t[fc] = {Vt[k])^^ with Vt[k] := exp 
where Vt[k] satisfies the SDE 

dVt[k] = Vt[k] (ik{t) dX{t) - -k{tf dt 



(157) 



(158) 



Follo wing (117) we again define a "density matrix" Xt[k]j and use Ito calculus 
(118) with the above equation for Vt[k] and the stochastic density matrix equa- 
tion for diffusive processes (156) to find that Xt[k] obeys an equation and initial 
conditions identical to ( |143 ). Thus the characteristic functionals <i>f[fc] defined 
in Eqs. (141) and ( |157D , respectively, are identical. 

As an example, we readily derive the mean homodyne current 



(159) 



/h°-(t)) ^ (/'--(i)),, = -^^^$^[fc]|,.o = T^sxpii) 



and the stationary homodyne correlation function whose Fourier transform gives 
the spectrum of squeezing H 



r°'^{ti)r°"'{t2) 



= H) 



Trs 



(/''°'°(il)/'^°'"(t2)) 



5k{t) Sk(t 



xe 



£(*2-ti) 



(cp + pc^)] +6{t2-h) 



(160) 



where the last term is the shot noise contribution. 

As a final comment we note that a completely different interpretation of 
the stochastic Schrodinger equations of the type (153) have been proposed in 



31 



Refs. [p4, pq, p^ in connection with dynamical theories of wave function re- 
duction. The assumption is that the reduction of the wavefunction associated 



with a measurement is a stochastic process and an equation of the type (154) 
is postulated. 

There have been numerous applications of the diffusive SSE. A beautiful 
discussion of squeezing in a parametric amplifier in terms of trajectories can be 
found in Carmichael's book M. Illustrations for the decay of Schrodinger cats 
have been given in |l^, |3|, |3^, Q. A theory of quantum feedback has been 
developed by Milburn and Wiseman ||l4| . 

0.5 Application and Illustrations 

During the last years we have seen numerous, mostly numerical applications of 

the SS E a iioi om, 111, m, [HI §1 0, §,11,1111 n, g, n. 



BOL p]J, p3, p3| . Roughly speaking, these calculations can be divided into illus- 
trations where single trajectories illustrate certain features of the physics one 
wishes to discuss, and as a numerical simulation method for situations where a 
direct solution of the master equation is not feasible. Below we discuss a few 
examples from our recent work. Our first example is state preparation in ion 
traps through quantum jumps; the second example is a discussion of simulations 
in laser cooling, and we conclude with remarks on possible application to the 
problem of decoherence in quantum computing. 

0.5.1 State Preparation by Observation of Quantum Jumps 
in an Ion Trap 

Ions traps provide a tool to store and observe laser cooled single ions for an 
essentially unlimited time. Experiments with single trapped ions thus represent 
a realization of continuous observation of a single quantum systemm. the context 
of quantum optics. For a review on ion traps and application in quantum optics 



we refer to |44, 4 



Quantum Jumps in Three-Level Atoms 

The probably best-known example is the problem of "observation of quantum 
jumps in three level system" [^ (for a review and references see [Q). The 
system of interest involves a double-resonance scheme where two excited states 
|e) and \r) are connected to a common lower level \g) via a strong and weak 
transition, respectively. The fluorescent photons from the strong transition are 
observed. However, an excitation of the weak transition where the electron is 
temporarily shelved in the metastable state \r) will cause the fluorescence from 
the strong transition to be turned off. It is, therefore, possible to monitor the 
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quantum jumps of the weak transition via emission windows in the (macro- 
scopic) signal provided by the fluorescence of the strong transition. Experi- 
mental observation of this effect has been reported in Refs. |El|, and various 
theoretical treatments of this effect have been published. From a theoretical 
point of view, the problem is to study the photon statistics of the photons emit- 
ted on the strong line of the three-level atom. A treatment based on continuous 
measurement theory was given by Zoller et al. [Q and Barchielli [^ (see also 
[ [46| ) : in Ref . |Q the probability density for the "emission of the next photon 
on the strong line" (delay function) was calculated, and for the first time a sim- 
ulation of the a posteriori dynamics was described to illustrate the conditional 
dynamics of a single continuously monitored quantum system. 
The master equation for a three-level atom has the form 

p=-i{H,sp-pHl^)+J,p + J^p i=Cp) (161) 

where jTsP :— \g){g\ ^sPee and jTwP '■— \g){g\ ^wPrr are the "recycling operators" 
for the atomic electron on the strong and weak transition, respectively. The 
probability density for the emission of a photon on the strong line at time t, 
given the previous photon on the strong transition was emitted at tr is, according 



to Eq. (106) 



c(s,i|s,t„) = pll'^\s,t\p,{t„))/dt 

Trs{Jse^^^^'^^'''-^JsP} 
TrsiJsp} 



(162) 



where Pc{tn) = \g){g\ is the density matrix at time i„ after emission of a photon 
on the strong line at time i„. We note that this conditional state is independent 
of the previous history of photon emissions; this is due to the fact that a quan- 
tum jump always prepares the atom in the ground state \g). In addition, we 
have summed over an arbitrary number of possible transitions on the weak line 



which replaces S{t, i„) -^ exp (£ — Js){'t — t„) in Eq. (102), p is the steady state 



density matrix. It is obvious from ( |l6l| )[p8| Eq. (162) is a sum of exponentials 
with different decay constants, reflecting the different time scales for the exci- 
tation and decay on the strong and metastable transition, which - when these 
distributions are simulated - manifest themselves in the periods of brightness 
and darkness in the fluorescence on the strong transition. In particular, the 
density matrix conditional to observing an emission window on the strong line 



when the last s-photon was emitted at time t„ is from Eq. (101) 



p,(t)-e(^--7=)(*-*'-)p/Trs{...}^|r)(r| (r,(t - t„) » 1) , (163) 

i.e. observation of a window in a single trajectory of counts corresponds to a 
preparation of the electron in the metastable state \r) ("shelving of the electron"). 
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Fock states in the Jaynes-Cummings model 

The principle of state preparation by quantum jumps can be extended to more 
complicated configurations. As a second example, we discuss the preparation 
of Fock states in a Jaynes-Cummings Model (J CM) based on the observation of 
quantum jumps |Q (see also Q). The JCM describes a harmonic oscillator 
strongly coupled to a single two-level atomic transition. The Hamiltonian is 

H — hujfa^a + —huJeg<Jz + hg{<j^a + aV_), (164) 

where a' and a are creation and annihilation operators for a mode of the radi- 
ation field with frequency w^, and (j±^z ai'e the Pauli spin matrices describing a 



two-level atom with transition frequency ojeg- The last term in (164) describes 
the coupling of the field mode to the atom with coupling strength g. Dissi- 
pation can be included in this model by coupling the field mode and atom to 
independent heatbaths and using a master equation formulation, in which one 
introduces damping rates k and F for the field mode and atom respectively. Of 
particular interest is the strong-coupling regime in which g > k when quan- 
tum effects in the coupled oscillator-spin system are most pronounced. Here 
the spectroscopy of the system is best described in terms of transitions between 
the "dressed states," or eigenstates, of the Hamiltonian H. The ground state 
is \0,g) — \0)\g) and, on resonance (wo = ^f), the excited dressed states are 
|n, ±) = {\n — l)|e) ± \n)\g))/\/2 {n — 1,2, . . .), where \g) and |e) are the bare 
atomic ground and excited states respectively, and \n) are Fock states of the field 
mode with excitation number n = 0, 1, 2, . . .. The eigenenergies corresponding 
to the excited states are En± — fi[ujQ{n — 1/2) ± g\/n\, which show AC-Stark 
splitting proportional to y^. The two lowest transitions |1, ±) — > |0, g) give rise 
to a doublet structure, the "vacuum" Rabi splitting M, B3. 

Experimental realization of the JCM with dissipation has been demonstrated 
in the field of cavity QED, both in the microwave and optical domain [Q [l9[ . 

An alternative realization of the JCM is a single trapped ion constrained to 
move in a harmonic oscillator potential and undergoing laser cooling at the node 
of a standing light wave. Under conditions in which the vibrational amplitude 
of the ion is much less than the wavelength of the light (Lamb-Dicke limit), this 
problem is mathematically equivalent to the JCM with negligible damping of 
the oscillator (i.e. k = 0). In this configuration, preparation of the Fock state 
corresponds to preparation of a non-classical state of motion of the trapped ion 
with fixed energy nhuj, whereas in cavity QED the Fock state corresponds to 
a non-classical state of light with no intensity fluctuations and undetermined 
phase. 

As discussed in detail in Rcfs. |4^ the master equation for single two-level 
ion {|g),|e)} trapped in a harmonic potential and located at the node of a 
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standing light wave has the form 

d 
di 



d 



iya^a+ (-^)<^z - 2"''('^+ + <^-){a + a''), p 



r 

+ — (2cr_pcr+ — 0"+cr_p — pcr+CT-). (165) 

where the oscillator now refers to the ion motion with frequency v (phonons), 
A is the detuning of the laser from atomic resonance, and the atom - oscil- 
lator coupling is proportional to the Rabi frequency Q times the Lamb-Dicke 
parameter 77. Eq. (|l65| ) is valid to lowest order in rj — nao/X <^ 1 (Lamb- 
Dicke limit) with Uq the size of the ground state wave function of the trap and 
A the wavelength of the light. Making the associations —A ^^ wq, 1/ <-* ojf, 
and g <-> ij^/2, we have a clear analogy with the damped JCM. The additional 
terms proportional to a-i-a* and cr_a are usually dropped in the optical regime 
on the basis of the RWA. In our instance, such an approximation requires that 
f, A ^ ri, r, \v— A|, which is consistent with the sideband cooling limit implied 
above. A significant feature of the JCM produced by this configuration is that 
dissipation in the system is entirely due to damping of the two-level transition, 
while the oscillator is undamped (i.e. k = 0). Further, the effective coupling 
constant rj^/2 depends on the Rabi frequency and, in contrast to CQED, can 
be adjusted experimentally to satisfy the strong coupling condition T < 77J7/2. 

The spectroscopic level scheme of the Jaynes-Cummings ladder of the strongly 
coupled ion + trap system is shown in Fig. ^. One possible approach to the 
observation of this level structure is a measurement of the probe field absorp- 
tion spectrum to a third atomic level \r), very weakly coupled to the otherwise 
strongly coupled |5)-|e) transition, analogous to level schemes used for the ob- 
servation of quantum jumps discussed above. We note that as a consequence 
of the unequal spacing of the energy levels in the JCM, the probe laser is only 
resonant with the transition frequency between a single pair of levels, and thus 
will only excite the system to a single state |n, r). 

This ability to selectively excite a particular transition, together with the 
state reduction associated with the observation of quantum jumps, offers the 
intriguing possibility of generating number states of the quantized trap motion. 
This follows from the fact that the probe laser exciting transitions to the states 
\n, r) interacts only with the atomic ground state contribution to the partic- 
ular dressed state |n, ±) being excited. Given that we are able to distinguish 
spectroscopically between the different maxima characterizing the absorption 
spectrum (so that we can identify the dressed state being excited), observation 
of an emission window in the presence of a weak coupling to the state \n,r) 
will tell us with certainty that the vibrational state of the ion is \n) (i.e. this is 
the vibrational state occurring with \g) in the dressed state |n,±)) and that we 
have produced a Fock state of the quantized trap motion. 

An obvious consequence of having the freedom to choose which transition 
is excited is the ability to choose the Fock-state that is to be produced. As 
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Figure 3: Level scheme of ion-trap systera. (a) Bare states \n,g),\n,e): The 
laser-atom detuning is lol — ujq = —v. Arrows pointing to the left indicate the 
resonant laser coupling with strength rfl\ arrows pointing to the right corre- 
spond to non-rotating-wave terms. There are no \n^g) — |n, e) couplings since 
the ion is at the node of the standing wave, (b) Dressed states |n,±) resulting 
from the laser coupling. B^ 

an example, Fig. |4|(a) shows the simulation of an experiment with a trapped 
three-level ion and the calculated random telegraph signal due to the probe 
excitation. Shelving of the electron in the state |2, r) are indicated when emission 
windows appear in the fluorescence of the strongly coupled atom-trap system. 
Thus a Fock state of the trap motion with n = 2 is prepared during these 
dark periods. Fig. ^a) shows the number of photons as observed in a real 
experimental situation where the integration time constant is long compared 
with the decay time of the strongly coupled system. To highlight the internal 
dynamics, we evaluate the temporal evolution of the mean quantum number (n) 
and the corresponding entropy of the system for a different set of parameters 
and over a timescale in which system relaxation is clearly visible. This is shown 
in Fig. H(b) and Fig. ffl(c) which demonstrate that after a quantum jump to 
the state |r) the system's entropy is suddenly reduced to zero (indicating the 
Fock state). After the ion returns to its strongly coupled states |n, ±) the mean 
value (n) approaches its thermal equilibrium value (here set by the mean photon 
number N characterizing the broadband thermal noise field) and the entropy 
increases to its steady state value. 
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0.5.2 Wave Function Simulations of Laser Cooling: Ap- 
plications to Quantized Optical Molasses 

One of the prime examples of wave function simulations in recent years is 
application to laser cooling of neutral atoms (for a review see JS^])- Exam- 
ples are description of quantized optical molasses in one-, two- and three- 
dimensional configurations and fluorescence and probe transmission spectra 
[ p^ p^ , 24, ^, |2^, ^, |2|, ^. Here we give a brief summary of our work 



on the spectrum of resonance fluorescence from quantized optical molasses, in 
particular the comparison between theory and experiment on the fluorescence 
spectrum as measured by Jessen et al. |5^ (see also [|2[ |3|), and more recent 
work on a fully quantum mechanical description of diffusion in optical molasses 
employing wave function simulations |B3 . 

Quantized Atomic Motion in Optical Molasses 

Laser cooling is typically accomplished in optical molasses, employing a config- 
uration of counterpropagating laser beams. The role of the laser is twofold: it 
provides a damping mechanism, and leads to the formation of optical poten- 
tials (the ac Stark shift of the atomic ground states) pM. The physical picture 
underlying laser cooling and spectroscopy of one-dimensional (ID) molasses is 
illustrated in Figs. 0, H and M for an angular momentum Jg = 1/2 to Jg = 3/2 
transition. Figure Mshows the atomic configuration for this model pR p& . We 
consider two counterpropagating linearly polarized laser beams with orthogo- 
nal polarizations, so that the positive frequency part of the electric field has a 
position dependent polarization 

e|^+^ (z, t) = £{sm{kz)e+i + cos(fcz)e_i)e-''^* (166) 

where lo is the frequency, k = 2tt/X the wave vector with A the wavelength 
of the laser light, and e±i are spherical unit vectors. For red laser detun- 
ings A = Lu ~ LUeg and low laser intensities, i.e. small saturation parameter 
s — rj^/2/(A^ + r^/4) << 1 with fl the Rabi frequency and F the spontaneous 
decay width, the Stark shifts of the two \g-) and \g+) ground states will form 
an alternating pattern of optical bipotentials U±{z). This is shown in the up- 
per part of Fig. 0: due to the large Clebsch-Gordan coefficients for the outer 
transitions (see Fig. pf), minima will occur for the state \g±) at positions with 
pure cr -light. In addition, spontaneous emission causes transitions between 
these potentials via optical pumping processes. In the semiclassical picture of 
Sisyphus-cooling [p6[ one considers an atom moving on one of these potential 
curves, say U-{z). Transitions to the other potential C/+(z) then occur prefer- 
entially from the tops of U^{z) down to the valleys of U+{z), so that on the 
average the atomic motion is damped. Quantum mechanically, laser cooling can 
be understood as optical pumping between the quantized energy levels (band 
structure) in the limit of level separations ?iWosc much larger than the optical 
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Figure 4: Simulation of quantum jumps as a function of time. The probe laser 
is tuned to |2, — ) -^ |2,7-). (a) Fluorescence intensity on the |e) - \g) transition 
(fiw = 2 X lO^^Fs, Fw = 10^'^Fs). (b) Time evolution of (n), and (c) the 
entropy, S = -Tr{pln(p)}, of the system {fl^ = O.lFs) Q- 
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pumping rate 70 = sT (which imphes large laser detunings) [ p5| . In the time 
domain this condition corresponds to a situation in which an atomic ccntcr-of- 
mass wave packet in the potential undergoes many oscillations with frequency 
Wosc before an optical pumping process occurs pa . An example of a band struc- 
ture is shown in the upper part of Fig. 1^. As a result of laser cooling the atom 
will occupy the lowest energy levels, and will thus be strongly localized. The 
lower part of [7| shows the corresponding localization of the atom in minima of 
the U±{z) potentials. Transitions between the vibrational states will manifest 
themselves as sidebands (Raman transitions due to optical pumping) in the 
atomic absorption and emission spectra [^ |2| p3[ . 

The basis of a theoretical discussion of laser cooling is the solution of the 
Generalized Optical Bloch Equations for the atomic density matrix comprising 
both the internal and external (center-of-mass) degrees of freedom |54|. The 
corresponding master equation for the density matrix p of the ground state 
manifold and ID motion in the z-direction is |23, pa] {h — 1) 



P= - T'{Hc«{z)p - pHcs{z)^) 

n + l 



J2 duN„{u)B„{z)e-'^''^pB„{z)^e'''''^. (167) 

/T -^ — 1 



N„(u) is the angular distribution of spontaneous photons emitted in the di- 
rection u — cos{ks, 6z) with polarization a = 0,±1, and 70 = sT/2 is the 
photon scattering rate. The first two terms on the right-hand side of the master 
equation involve the non-Hermitian atomic Hamiltonian 

Hosiz) = ^-iUo + ^70) E BAz^B^iz) (168) 

describing the motion of the atomic wavepacket with kinetic energy p^/2Af in a 
multicomponent optical potential with depths determined by Uq = s|A|/2. The 



real part of the potential in (168) gives rise to quantized energy levels (band 
structure), while the imaginary part describes a loss rate due to optical pumping. 
The operators B„{z) correspond to Raman transitions between the ground state 
levels by absorption of a laser photon and subsequent emission of a spontaneous 



photon with polarization a, see Ref. [g5||. The master equation (167) is of the 



Lindblad form with an infinite number of channels (the u integration over the 



z-projection of the emitted photon wave vector in Eq. (167)) 

We are interested in the spectrum of resonance fluorescence, emitted along 
the z-axis, with frequency lj' and polarization a. This spectrum is proportional 
to the Fourier transform of the stationary atomic dipole correlation function 

cAt - to) = {BAtV e'"'-'^'^ B^to) e-''^^^^*")) . (169) 
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According to Section 0.4.4 and Ref. 123] we have 



c,(i-to)=Trs{i?ie''=^^p(+)(i)} (i > io) ; (170) 

where the first order perturbed density operator p^^' (t) obeys the same master 
equation as the density matrix but with a different initial condition 

p(+)(io)-i?.(5)e-*'=^*p(to). (171) 

Physically, this simulation procedure can be interpreted as the following 



computer experiment ||13|, |2C|, g5| (for a more formal derivation see Sec. 0.4.4). 
We surround the atom by an (infinite number of) unit efficiency photodetectors 
covering the Att solid angle. The simulated "click" of the photodetectors deter- 
mines the time, polarization and direction of the emitted photon (see Fig. o a). 
This is the simulation of the density matrix of the laser cooled atom. To simu- 
late the fluorescence spectrum, we place a Fabry-Perot in front of one of these 
photodetectors and "measure" the transmitted photon flux. This gives the spec- 
trum as a function of the tuning of the resonance frequency of the Fabry-Perot 
transmission (see Fig. H b) . 

In the experiment by Jessen et al. the resonance fluorescence spectrum from 
Rb^^ atoms in ID molasses is observed. The Rb atoms are driven on the 55*1/2 
Fg = 3 -^ 5P3/2 i^e = 4 transition in a lin _L lin laser co nfigu ration described 



above. In this case a direc t sol ution of the master equation (|167| ) to calculate the 



autocorrelation function (17C) is impractical due to the large dimensionality of 
the density matrix equation which involves N'^ elements (N = N^x x Nj^t where 
Nint is the number of internal, and N^x the number of (discretized) external 
degrees of freedom). For ^^Rb we have N — 448 on a Fourier grid with 64 points 
corresponding to momenta up to zL32hk [E5|; 10000 wavefunction realizations 



are needed for convergence. Adapting the formalism of Sec. 3.4.5 a simulation 



of the master equation (167) consists of propagation of an atomic wave function 
ipcit) with the non-Hermitian (damped) atomic Hamiltonian (|l6q ) interrupted 
at random times by wave function collapses 

i,,it + dt) = e''''^'B,iz)Mt), (172) 

and subsequent wave function renormalization. The Schrodinger equation for 
tpcit) describes the time evolution of the atomic wave packet in the periodic 
optical potential, and its coupling to the laser driven internal atomic dynamics. 
The times of the "quantum jumps" are selected according to the delay function 



ciu',a,t) = llv/70 N^u') B4z)Mt)f (173) 

which gives the probability for emitting a spontaneous photon at time t, with 
momentum k'^ = ku' along the z-axis and polarization a. The "quantum jump" 



(172) corresponds to an optical pumping process between the atomic ground 
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Figure 5: (a) The "experiment" simulated in the Monte Carlo wave function 
description of optical molasses: the atom is surrounded by photo detectors cor- 
responding to an angle resolved detection of the emitted photon h ks- (b) 
Simulation of the fluorescence spectrum of the emitted light: the atom is "ob- 
served" with photodetectors and the spectrum is obtained by filtering one output 
channel with a Fabry-Perot interferometer p5| |. 

states, including the associated momentum transfer to the atom. Averaging over 



these wave function realizations gives the density matrix according to Eq. (121). 



The dipole correlation function (|170|) can be simulated following the ap- 



proach outlined in section ].4.4 [llSL ECj, 25 . The perturbed density matrix 



p(+)(f) in Eq. (|170|) can be interpreted as a first order response to a "delta- 



kick" at time t = to, represented by the initial condition (171,132|). A sim 



ulation is obtained by introducing a "perturbed" wave function ipc (t) which 
obeys the Schrodinger equation for tpcit) but now with initial condition [compare 
Eq. dmlp^ )] 

i^^+Ht = to) - e-^'=^^i?,(z)^e(io), (174) 

and quantum jumps of ipc (t) dictated by the wave function ipcit) according to 



the delay function (173). The dipole correlation function is 



cAt-to) - {mt)\B,{z)^e<'\ipi+Ht))/\mt)r)^^ (175) 

where the angular brackets indicate averaging over both quantum jumps and 
initial times to- 

The periodicity of the atomic Hamiltonian ( |l6q ) in space, and assuming an 
infinite periodic molasses, allows one to propagate the atomic wave packets as 
time dependent Block functions 

Mz,t)^^e''''ug{z,t) (176) 

V 27r 

with q £ [—k,k] a quasi momentum in the first Brillouin zone, and Uq{z,t) = 
Uq{z + A/2,i) a periodic multicomponent atomic wave function. The Hamil- 
tonian evolution due to Hcs preserves q, while quantum jumps cause changes 
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lg-l/2> lgl/2> 



Figure 6: Atomic level scheme and Clebsch-Gordan coefficients for a Jg = 1/2 
to Je = 3/2 transition in a lin _L lin configuration [E5| . 



between families of Bloch functions, q -^ q' . In practice, we propagate the Bloch 
function Uq{z,t) on the unit cell of the lattice z e [0, A) using a split-operator 
Fast-Fourier Transform method |23] . 

Figure ra compares the resonance fluorescence spectrum for cr+ polarized 
light obtained by simulation (solid line) with the experimental data of Jessen 
et al. pu (crosses) for laser intensities, detunings etc. taken directly from the 
experiment with no adjustable parameter. The theoretical spectrum of Fig. 
was obtained by convolving the ab initio spectrum with the experimental resolu- 
tion. The central line is scattering at the laser frequency while the first red and 
blue sidebands correspond to Raman transitions between adjacent vibrational 
bands in the optical potential. The asymmetry of the red and blue sideband 
intensities reflects the populations of the vibrational levels, and from the excel- 
lent agreement between theory and experiment we infer that the wave function 
simulation reproduces the experimental temperature of the atoms [plf . Laser 
cooling accumulates atoms predominantly in the lowest vibrational states and 
in the Mg = ±3 potentials. This spatial localization of atoms - on a scale small 
compared with the laser wavelength - suppresses optical pumping transitions 
between different vibrational levels n ^ n\ which is responsible for a narrowing 
of the lines in the optical spectrum (Lamb-Dicke narrowing). A broadening 
mechanism is present for the sidebands, due to the anharmonicity of the opti- 
cal potential. This leads to different transition frequencies for n — ^ n ± 1 (by 
approx. Er/K). For the present parameters this anharmonicity is not resolved 
even in the unconvolved spectrum. 
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Figure 7: In the upper panel the optical bipotentials and bandstructure of 
the atom are plotted as a function of position z for a 1/2 to 3/2 transition. 
The potential depth is C/q = lOOE'ij. The excited states were adiabatically 
eliminated. We schematically indicate the two Raman processes between the 
ground and excited states which lead to the red and blue sidebands in resonance 
fluorescence. In the lower panel we show the spatial distribution of the atoms 
in the |g+)-state (solid line) and |5_)-state (dashed line) for U^ = 100i?_R and 
7o = ^/^Er. The atoms are localized in the valleys of the corresponding optical 
potentials |25[. 
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Figure 8: Spectrum of resonance fluorescence as a function of the frequency 
u. The solid line is the theoretical spectrum [g5| convolved with a Lorentzian 
corresponding to a finite detector width of 3.8 kHz and a Gaussian with width 
20 kHz (residual Doppler broadening). Crosses are the experimental result of 
Jessen et al. [^ . The inset shows the total spectrum. For parameters we refer 
toRef. [El. 



Localization by Spontaneous Emission 



As outlined in the context of Eq. (122) there is no unique way of decomposing a 
given master equation ( p8| ) to form quantum trajectories ipc{t): since the reser- 
voir measurement may be performed in any basis. This statement is equivalent 
to noting that Eq. (| 



g) is form invariant under the substitution (122). Different 
sets of jump operators {cj} not only lead to a different physical interpretation of 
trajectories, but an appropriate choice of Cj may be crucial for the formulation 
of an efficient simulation method for estimating the ensemble distribution [Q . 
We will illustrate these ideas in the context of the q uanti zed motion of an atom 
moving in optical molasses for the master equation (167). 

Simulation methods for modeling spontaneous emission according to the 
master equation (167) have assumed an angle resolved detection of the photon 
as illustrated in Fig. 1(a). For a one dimensional sy stem with adiabatically 
eliminated excited state this gave the jump operators ( |172[ ), 



v/7o^.(w)e-^'="^B,(z) (-l<ii<+l,(7-G,±l) 



(177) 



In a simulation each quantum jump gives information on the direction of the 
emitted photon. Alternatively, we could observe the fluorescence through a lens 
(Fig. ^(b)). The wave function simulation in this case is equivalent to the direct 
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simulation of a Heisenberg microscope [^: it is not possible to distinguish 
different paths ui and M2 (Fig. ^(b)) which may be taken by the photon so the 
decay operator does not generate a unique recoil. Instead, information about 
the position coordinate is provided by each emission event, and as a result the 
wavefunction is l ocali zed in position space. Applying a Fourier transform to the 



operators in Eq. (177) to model the action of the ideal lens gives the new decay 



operators for one dimension 



du v/7oA^.(w)S,(i)e-''="(*-''^/2) (^,^o,±l,...). (178) 



-1 



For angle resolved detection, u labels a continuous but bounded set of operators, 
so that in the conjugate basis, v can be any integer and indexes an infinite set 
of operators at discrete points. The integral can be evaluated for c^^q „ to give 
a localized function centered at the origin and the rest generated by t ransl ation 



by m ultiples of ±A/2. To prove that both sets of jumps operators (177) and 



( 178 ) give rise to the same a p riori dynamics, we consider the following identity 



for the recycling term in Eq. ( 167 ) 



1 



J ,poc Jdu (TiV^B^e'""^) p ( Viv^B^e''^"^) ^ = 

-1 

= f2 ( [du y^^B,e^fe-(-'^V2) ]pl fdu' y^Mfl5^e./cn'(£-.A/2) 

where we have used 

- 00 
S{u - ^i') = 2 J2 e*'^'''"-"'' (-1 <u,u' < +1) . (180) 

z^— — 00 

As an application of the new simulation method, we consider a fully quantum 
mechanical treatment of atomic diffusion in optical molasses for the lin _L lin 
laser configuration described above |42| . The calculation of quantum diffusion 
using the angle resolved detection approach is difficult because the wave func- 
tion spreads out during the coherent propagation. In contrast to a description 
of an infinitely extended periodic molasses pa, we have here an intrinsically 



non-periodic problem. Applying the localizing jump operators (178) allows us 
to make use of a greatly reduced basis set which we allocate dynamically to 
follow the atom. One representative trajectory illustrating the random walks 
is shown in Fig. |l^(a) where we plot the expectation value of the spatial co- 
ordinate as a function of time. The long periods when the position does not 
change appreciably correspond to sub-barrier motion when the total energy of 
the atom is below the threshold given by the maximum of the optical poten- 
tial. Energy fluctuations allow eventually the atom to overcome the potential 
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atom photodetector jj^ 

atom lens photodetector 



Figure 9: Schematic illustration of the two different measurement bases for the 
detection of spontaneous emission [Q . 

barrier as indicated by the dashed line in Fig. M(b). It may then travel over 
several wavelengths until it is trapped again. For a more complete discussion 
of the spatial diffusion coefficient we refer to Ref. [Q. Other candidates for 
application of this method are cold collisions Q . 

0.5.3 Quantum Computing, Quantum Noise and Contin- 
uous Observation 

Quantum computers (QCs) (for a review see |59[) promise exponential speedup 
for certain classes of computational problems [|60| at the expense of exponen- 
tial sensitivity to noise [pll p3, ra, pM. Thus the analysis and suppression of 
quantum noise is a fundamental aspect of any discussion of the practicality of 
quantum computing. The theory of quantum noise and continuous measure- 
ment discussed in the preceding sections provides the theoretical basis for such 
an analysis. 

Computation is a process which transforms an initial state to a final state by 
the application of certain rules. These input and output states are represented 
by physical objects, and computation is the physical process which produces the 
final from the initial state. Quantum computers are physical devices which obey 
the laws of quantum mechanics, in the sense that the states of the computer 
are state vectors in a Hilbcrt space, and the processes are the unitary dynamics 
generated by the Schrodinger equation on this Hilbert space. The distinctive 
feature of QCs is that they can follow a superposition of computational paths 
simultaneously and produce a final state depending on the interference of these 
paths |p9| . Recent results in quantum complexity theory, and the development 
of some algorithms indicate that quantum computers can solve certain problems 
efficiently which are considered intractable on classical Turing machines. The 
most striking example is the problem of factorization of large composite numbers 
into prime factors iQ, a problem which is the basis of the security of many 
classical key cryptosystems: to factor a number N oi L digits on a classical 
computer requires a time approximately proportional to exp L^'^ for the best 
known algorithms, while Shor has shown that for a QC the execution time grows 
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Figure 10: (a) Expectation value of the spatial coordinate for a single trajectory 
versus time, (b) Corresponding (kinetic plus potential) energy expectation value. 
Note the coincidences between above barrier energies and long flight periods over 
many wavelengths Q . 
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only as a polynomial function L? pO| . 

Realization of the elements of a quantum computer 

The basic elements of the QC are the quantum bits or qubits [^ . In a classical 
computer a bit can be in a state or 1, while quantum mechanically a qubit can 
be in a superposition state It/;) = co|0) + ci|l) with |0) and |1) two orthogonal 
basis states in a two-dimensional state space 7^2 ■ Quantum registers are defined 
as product states of L qubits, 

\x) = \xl-i)l-i ■ ■ ■ No)o e ^2 ® . . . ^2 (181) 

with and a;„ — 0, 1, and x — X]n=o ^n^" the binary decomposition of the number 
X. There are 2^ different states and the most general state of a quantum register 
is the superposition 

£={0,1}-^ 



In a QC the physical objects are state vectors * in the product space (182), 
and the computation is the Schrodinger time evolution V'in — > V'out = Uipi-a with 
U a unitary time evolution operator. Note that unitarity of U implies that the 
computation is reversible. 

Physical problems in realizing a QC are: 

(i) A physical realization for the qubits must be found. Examples are internal 
states of atoms, or photon number states or polarization states of photons. 

(ii) We must be able to erase the content of the register and to prepare the 
initial superposition state i/'in- 

(iii) The unitary evolution operator U must be implemented. The unitary trans- 
formation can be decomposed into a sequence of steps involving the condi- 
tional dynamics of a few qubits (quantum gates) |5^ . 

In particular, it has been shown that any operation can be decomposed into 
controlled-NOT gates between two qubits and rotations on a single qubit, 
where a controlled-NOT is defined by 

Ci2:|ei)|e2)^|ei)|eiee2) £1,2 = 0,1 (183) 

where ® denotes addition modulo 2. The question then is to find physical 
processes (i.e. a Hamiltonian) which implement these quantum gates. 

(iv) To read the state of the register ■0out at the end of the computation we 
must implement and perform a von Neumann state measurement. These 
measurements are destructive, i.e. irreversible. 

(v) In practice the central obstacle to building a QC is the fragility of macro- 
scopic superpositions with respect to decoherence due to coupling of the 
QC to the environment, and the inaccuracy inherent to state measurements 
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[nil p3, p3l pjl . Schemes and experimental realizations must be found which 
minimize these sources of decoherence. 



Realization of quantum gates has been discussed in Refs. |^, |65| |66|, |6j, |68| . 
Cirac and ZoUer |g^ have proposed the implementation of a QC using laser 
cooled ions in a linear ion trap. The qubits are stored on a metastable atomic 
transition, while the ions undergo collective quantized motion in the linear trap 
[ [70| , and interact individually with laser light. Exchange of phonons stimulated 
by laser pulses allows the implementation of quantum gates. The distinctive 
features of this system are: it allows the realization of quantum gates between 
any set of (not necessarily neighboring) ions; in comparison with other systems, 
decoherence is small, and the final readout of the output state can be performed 
with essentially unit efhciency employing the methods of quantum jumps. A 
first experimental implementation for quantum gates along these lines has been 
reported from the NIST Boulder group [[71[ . A QC based on optical CQED 
has been proposed by Pellizzari et al. ]72| . It is based on a string of atoms 
trapped inside a high Q optical cavity. Long lived Zeeman coherences of ground 
state atoms represent the qubits. Quantum gates between atoms are realized 
via exchange of a photon through the single quantized cavity mode W^ . As in 
the case of ion traps, quantum gates can be performed between any two qubits. 
Sources of incoherence in this model are cavity decay and spontaneous emission 
of atoms from the excited state during the gate operation. However, spontaneous 
emission can be significantly avoided by performing the gate operation between 
two atoms as an adiabatic passage process via a dark state of the strongly 
coupled atom+cavity system |73|| ; and cavity decay is minimized by having a 
photon present only during the gate operation. 

Errors and their correction 

There are two sources of errors in quantum computing. First, there are sys- 
tematic errors, for example, due to uncertainties in adjustment of parame- 
ters. Second, there is decoherence as a result of coupling to the environment 
[ pi] |6^ , ^, |64|. The general state of a register of a QC is described by a density 
matrix 

2^-1 2^-1 

/'W = E E \^)p^,^'it)i^'\ ■ (184) 

x,=0 x'=0 

Thus as a result of the contact with the environment there will be changes in 
populations (the diagonal density matrix elements), and decay of the coherences 
(the off-diagonal density elements) on a time scale determined by the coupling of 
the QC to the environment. The power of quantum computing results from in- 
terference between quantum paths, consequently a loss of coherence will destroy 
the quantum parallelism. Unruh [pl| has shown that the computer-environment 
coupling leads to an exponential error rate with the number of qubits L, and 
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it has been argued that this dependence will be the decisive factor to limit the 
size of problems which can be solved on the QC. What these limits are will 
be determined by fundamental limits of specific theoretical proposals for the 
decay of qubits and noise in quantum gates, and by technological and experi- 
mental progress. Specific models for dissipation have been studied in references 

[imiiii. 

The above discussion of decoherence assumes that the environment is traced 
out, i.e. in the language of continuous measurement theory "the results of the 
measurements in the out channels are not read." This is the situation where the 



time evolution of the density matrix (184) is described by the master equation 



(a priori dynamics). It is of interest to study the evolution of a QC which 
is continuously monitored to detect the decay (quantum jump), and thus an 



incoherent step in the evolution (compare nonlinear evolution Eq. (103, and 
the quantum jump described by Eq. (|l05| ) |7^ . After detection of a decay the 
computer can be restarted; on the other hand, by restricting the ensemble to 
those computer runs with no decay, we improve the statistics for the read out. 
In the CQED example given above this could be a cavity decay, or a photon 
from a spontaneous emission of one of the atoms representing the qubits. Note, 
however, that - even when no decay is detected in a specific run - information is 
gained by this "non-observation of a decay" and this leads to a distortion of the 
superposition of the qubits. To illustrate the effects of damping and continuous 
observation, we consider the following simple example. Let us assume that we 
store a superposition (qubit) |V'(t = 0)) = co|0) +ci|l) on an atomic transition 
with ground state |0) and excited state |1), where the upper state decays with 
a radiative rate T. At time t = we prepare the qubit in the superposition 
\^p(t = 0)) — co|0) + ci|l). If we do not monitor the decay of the qubit (a priori 
dynamics), the time evolution of the qubit is governed by the solution of the 
master equation: 

\^ ( i = 0))(V(i = 0)|^exp/:t|V^(0))(V'(0)| (185) 

= |conO)(0| (1 - e-^') + |cini)(l|e-" + (cic5|l)(0| + H.c.) e"*^* . 

We note that the initial pure state develops into a mixture, until the atom finally 
has decayed to the ground state. If we monitor the system continuously (with 
unit detection efficien cy) a nd observe no decay (a posteriori dynamics), the state 
evolves according to (|l03| ) and remains a pure state; however the superposition 
is distorted by the decay (see also Ref. jig] ): 

mt^O)) ^ |V'(0>-e-'^-*|^(0))/||e-^^-*|^(0))|| 

= (co|0)+ci|l>e-^i^*)/||...|| (186) 

-^ |0) iort^oo . 

This unwanted dynamics of the qubits corresponds to our increase of information 
which we gain by observing (measuring) the system: by not registering a decay 
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our knowledge increases that the system is in the ground state. The distortion 
of the qubit is a result of the asymmetry of the decay. When we flip the qubit 
at time t/2, |0) ^-> |1) , and again at time t, we symmetrize the decay, 

m = 0)) ^ (co|0)e-^i^* + ci|l)e-^"^/||...|| (187) 

= co|0)+ci|l) = |^(t = 0)). 

Thus, the information which state has decayed has been lost, the superposition 
remains unchanged, and there is no effect of decay (in the subensemble with no 
decay) . The idea of a symmetric decay, which factors out in the time evolution 



(103) and thus drops out upon normalization, is readily generalized to a con- 
struction of gates. This problem is studied in more detail in Ref . 173] . It is an 
interesting question to what extent schemes can be found where the state of the 
register can be reconstructed after a quantum jump 174]. 
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.1 An exercise in Ito calculus 

As an exercise in Ito calculus we will go in some detail through the steps of 



deriving Eq. (|ll2| ) from ( |110| ). To simplify notation we confine ourselves to the 



case of a single channel Nc = 1- We will show that if tpc{t) obeys 

d4(t) ^[Adt + {B- i)diV(i)]4(i) (188) 

with A and B operators, then the normalized ijjc{t) obeys 

dMt) = [{A -{A + A^)J dt + {B/^{B^B)^ - l)dN{t)]i:rit) (189) 



with (. . .)g defined in Eq. (|11 
Proof: From 

^c{t + dt) = [i + Adt + {B- l)dN{t)]Mt) (190) 

the change of the norm in the time step dt is 

Uc{t + dt)f - \mt)f [1 + {A + At)^ dt + {B^B - i)^ dN{t)\ (191) 

and 
mt+dt)\\-^ - mt)r^[\ + {A + A^)^dt + {{B^B-l))jN{t)]-^'^ 

= mt)\\-\^-\{A + A^),dt + {l/^{B^^-l)dN{t)]. 
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This follows from 



1 



f{l + adt + bdN) = f{l) + f'{l)adt + f'{l)bdN+^f"{l)b^dN^ + ... 
= f{l) + ril)adt+[f{l + b)-fil)]dN 

valid for sufficiently well-behaved functions /. Finally, we have 

ij;^{t+dt) = [i + A dt + {B - i)dN{t)]^Pcit) X (192) 

[1 - 1 (A + At)^ dt + {l/^{BtB)^ 1) dNit)] 



which gives ( |189D 
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